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^ ' Gravitational collapse is one of the most fruitful subjects in gravitational physics. It 

\ is well known that singularity formation is inevitable in complete gravitational collapse. It 

5— ( , was conjectured that such a singularity should be hidden by horizons if it is formed from 

generic initial data with physically reasonable matter fields. Many possible counterexamples 
' to this conjecture have been proposed over the past three decades, although none of them 

has proved to be sufflciently generic. In these examples, there appears a singularity that is 
. not hidden by horizons. This singularity is called a 'naked singularity.' The appearance of 

' a naked singularity represents the formation of an observable high-curvature, strong-gravity 

" ^' ' region. In this paper we review examples of naked singularity formation and recent progress 

, in research of observable physical processes - gravitational radiation and quantum particle 

■ creation - from a forming naked singularity. 
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§1. Introduction: Review of naked singularity formation 



1.1. Introduction 

According to the general theory of relativity, which is a classical gravity theory, 
after undergoing supernova explosion in the last stage of its evolution, a star with a 
mass dozens of times larger than the solar mass will contract without limit, due to 
its strong gravity, and form a "domain" called a "spacetime singularity." Spacetime 
singularity formation is a very general phenomenon, not only in the gravitational 
collapse of stars of very large but also in physical processes in which the general 
theory of relativity plays an important role. In fact, it was proved by Hawking and 
Penrose that the appearance of spacetime singularity is generic, i.e., spacetime 
singularities appear for any spacetime symmetry. However, the singularity theorems 
of Hawking and Penrose only prove the causally geodesic incompleteness of spacetime 
and say nothing about the detailed features of the singularities themselves. For 
example, we cannot obtain information about how the spacetime curvature and the 
energy density diverge in a spacetime singularity from these theorems. 

Spacetime singularities can be classified into two kinds, according to whether or 
not they can be observed. A spacetime singularity that can be observed is called a 
"naked singularity" , while a typical example of a spacetime singularity that cannot 
be observed is a black hole. A singularity is a boundary of spacetime. Hence, in 
order to obtain solutions of hyperbolic field equations for matter, gauge fields and 
spacetimes themselves in the causal future of a naked singularity, we need to impose 
boundary conditions on it. However, we do not yet know any law of physics that 
determines reasonable boundary conditions on singularities. Therefore, the existence 
of a naked singularity implies behavior that cannot be predicted with our present 
knowledge. 

Is such a naked singularity formed in our universe? With regard to this question, 
Penrose proposed the so-called cosmic censorship conjecture. This conjecture 
represents one of the most important unsolved problems in general relativity. Its 
truth is often assumed in the analysis of physical phenomena in strong gravitational 
fields. There are two versions of this conjecture. The weak conjecture states that 
all singularities in gravitational collapse are hidden within black holes. This con- 
jecture implies the future predictability of the spacetime outside the event horizon. 
The strong conjecture asserts that no singularity visible to any observer can exist. 
It states that all physically reasonable spacetimes are globally hyperbolic. Unfor- 
tunately, no one has succeeded in giving a mathematically rigorous and provable 
formulation of cither versions of the cosmic censorship conjecture. 

If naked singularities are formed frequently in our universe, then they have a very 
important meaning in the experimental study of physics in high-energy, high-density 
regimes. Solutions of theoretical models exhibiting naked singularity formation may 
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be nothing more than unreaHstic behavior of toy models and cannot be considered 
as providing proof of the existence of naked singularities in our universe. However, 
some such solutions are likely to become important footholds from which we can 
advance research of spacetimc singularities. In fact, the possibility of naked singu- 
larity formation in a large linear collider has recently been suggested in the context 
of the scenario of large extra dimensions. In this paper we review recent progress in 
theoretical research on naked singularity formation and physical processes contained 
therein. In the remaining of this section, we outline examples of naked singularities 
in the case of spherically symmetric gravitational collapse. We also outline the hoop 
conjecture and research on axisymmetric and cylindrical collapse. In §2 we review 
recent works on nonspherical perturbations of spherical collapse and the gravita- 
tional radiation from a forming naked singularity. In §3 wc review quantum particle 
creation from a forming naked singularity. We use units in which G = c = h = 1. 

1.2. Spherical dust collapse 

In order to examine the validity of the cosmic censorship conjecture, it is clear 
that we have to deal with dynamical spacetimes. In this context, we often assume 
some kind of symmetries for the spacetime, because the introduction of symme- 
try makes the analysis much easier. Among such symmetric spacetime, spherically 
symmetric spacetimes have been the most widely studied partly because for them 
the analysis becomes much simpler due to the absence of gravitational waves, and 
partly because we can expect that some class of realistic gravitational collapse can 
be treated as spherically symmetric with small deviations from it. Here we review 
spherically symmetric gravitational collapse and the appearance of naked singulari- 
ties in spherically symmetric spacetimes. 

1.2.1. Spherically symmetric spacetime 

Before restricting matter fields, we present the Einstein equation for a general 
spherically symmetric spacetime. In spherically symmetric spacetime, without loss 
of generality, the line element can be written in diagonal form as 

Here we adopt a comoving coordinate system, which is possible for matter fields of 
type I. (See Ref. 7) for classification of matter fields.) In this coordinate system, the 
stress-energy tensor T'^ that is the source of a spherically symmetric gravitational 
field must be of the form 
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(1-2) 



where p(t,r), U{t,r) and n{t,r) arc the energy density, radial stress and tangential 
stress, respectively. If we consider a perfect fluid, which is described by 



T^'^ = {p + P)ufu'' + Pgt"', 



(1-3) 
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then the stress is isotropic, i.e., 

E = n = p. (1-4) 

^From the Einstein equation and the equation of motion for the matter fields, 
we obtain 



m! = 




(1-5) 


m = 


-AttER^R, 


(1-6) 


R' = 


Ru' + R'X, 


(1-7) 


E' = 


-{p + Ey-2{E-n)^, 


(1-8) 


m = 


|(i-i?'V2^ + i?V2-), 


(1-9) 



where m = m{t, r) is the Misner-Sharp mass, and the prime and dot denote partial 
derivatives with respect to t and r, respectively. 

Here we stipulate the existence of an apparent horizon, which is defined as the 
outer boundary of a connected component of the trapped region. The important 
feature of the apparent horizon is that, if the spacetime is strongly asymptotically 
predictable and the null convergence condition holds, the presence of the apparent 
horizon implies the existence of an event horizon outside or coinciding with it. If 
the connected component of the trapped region has the structure of a manifold with 
boundaries, then the apparent horizon is an outer marginally trapped surface with 
vanishing expansion. (See Ref. 7) for the definitions and proofs.) Along a future- 
directed outgoing null geodesic, the relation 

is satisfied, where the upper and lower signs correspond to expanding and collapsing 
phases, respectively, and we assume R' > 0. Therefore, in the expanding phase, 
there is no apparent horizon. In the collapsing phase, on a hypersurface of constant 
the two-sphere R = 2m is an apparent horizon. The region R < 2m is trapped 
{dR/dt < 0), while the region R > 2m is not trapped {dR/dt > 0). 

Here wc should describe singularities that may appear in spherically symmetric 
collapse. A shell-crossing singularity is one characterized hy R' = and R > 0, while 
a shell-focusing singularity is one characterized by i? = 0. Also a central singularity 
is one characterized by r = 0, while a non-central singularity is one characterized by 
r > 0, where r = is chosen as the symmetric center. It is noted that, since 

R^' ^ = 8tt{p - E - 2n) (1-11) 
R^'''R^^ = 647r2(p2 + e'^ + 277^) (1-12) 

are satisfied, the divergence of the energy density or stress directly implies a scalar 
curvature singularity. As for strength of singularities, two conditions are often used; 
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one is the strong curvature condition proposed by Tipler, and the other is the 
hmiting focusing condition proposed by Krolak. The former condition is stronger 
than the latter. (See also Ref. 11).) These conditions are defined in terms of the 
speed of growth of the spacetime curvature along a geodesic that terminates at the 
singularity. A shell-crossing naked singularity is known to be weak in terms of both 
conditions. It is expected that a spacetime with weak singularity can be extended 
further in a distributional sense, although it is not known how this is possible in 
general situations. 

1.2.2. Dust collapse 

First we consider a dust fluid, which is defined as a pressureless fluid. This 
model has been most widely studied, partly because of the existence of an exact so- 
lution, which is called the Lemaitre-Tolman-Bondi (LTB) solution, " and partly 
because it provides a nontrivial interesting example of naked singularity formation. 
For this reason, we describe this model in some detail here. 

We restrict the matter content of the model to a dust fluid. Therefore 

E = 0, (1-13) 

n = 0. (1-14) 

Then, Eqs. (l-5)~(l-9) can be integrated as 

m=^, (M5) 
F' 

P= o_o2o/ ' (1-16) 



u = 0, (1-18) 
^' = / + ^, (1-19) 

where the arbitrary functions -F = F{r) and l + /=l + /(r)>0 are twice the 
conserved Misner-Sharp mass and the specific energy, respectively. In Eq. (1-18), 
we have used the rescaling freedom of the time coordinate. This means that a syn- 
chronous comoving coordinate system is possible. Equation (1-19) can be integrated 
to yield 



VF \ F 
where R^{r) and G{y) are defined as 



(1-20) 



i?°(r) = i?(0,r). 



(1-21) 
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for < y < 1 
for y = 

for y < 0, 



the quantity [Q(-R)];^o is defined as 



[Q(i?)]^o=Q(i?)-Q(i?°), 



(1-22) 



(1-23) 



(1-24) 



and the upper and lower signs in Eq. (1-20) correspond to expanding and collapsing 
phases, respectively. Hereafter, our main concern is with the collapsing phase. 

Assuming that R is initially a monotonically increasing function of r, and rescal- 
ing the radial coordinate r, we identify r with the circumferential radius R on the 
initial space-like hypersurface t = 0. Then, regularity of the center requires 



/(O) 
i?(t,0) 
F{r) 



= 0, 
= 0, 

< GO at r 



0. 



(1-25) 
(1-26) 

(1-27) 



The solution can be matched with the Schwarzschild spacetime at an arbitrary radius 
r = Tfe if we identify the Schwarzschild mass parameter with F(rh)/2. 

It is helpful for later use to write down the solution for marginally bound collapse, 
i.e., the case of /(r) = 0. This solution is 



R 



4 J 



) (ts - tf'\ 



where ts = ts{r) is given by 



tsir) 



.3/2 



Then, R can also be written as 



R = r\l- 




(1-28) 
(1-29) 

(1-30) 



The matching condition with the Schwarzschild spacetime yields the relation between 
the Schwarzschild time coordinate T and the synchronous comoving time coordinate 
t at r = ri, as 



T = t-2\/2Mi? + 2Mln 

VR - V2M ' 

where an integral constant has been absorbed through the redefinition of T. 

Again we go back to general LTB solutions. These solutions allow shell-crossing 
singularities which may be naked. Hereafter, we concentrate on shell-focusing 
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singularities. Equation (1-19) implies that every mass shell labelled by r that is 
initially collapsing inevitably results in shell-focusing singularity. It is easily found 
that the time at which the shell-focusing singularity appears, is(r), and the time of 
the apparent horizon, tAHir), are given by 

t,{r) = '-^G(-^), (1-32) 



^/F \ F 

tAHir) =ts{r)-FG{-f). (1-33) 

Therefore, a shell-focusing singularity that appears at r > is in the future of the 
apparent horizon. 

A non-central shell-focusing singularity is not naked. Indeed, suppose that a 
light ray emanates from a shell-focusing singularity at some ri > 0, which is given 
hy t = t{r). Then, by continuity there must exist an e > such that for ri < 
r < ri + € a light ray with positive expansion is later than the apparent horizon 
and earlier than the shell-focusing singularities, since the apparent horizon is earlier 
than the shell-focusing singularities everywhere but at the center. This implies that 
< R{t{r),r) < F{r) and dR/dt{t{r),r) > 0. By Eq. (1-10) these relations lead to 
a contradiction. We thus conclude that shell-focusing singularities (except possibly 
for central shell-focusing singularities) are not visible to an observer. Therefore it is 
sufficient to consider central shell-focusing singularities in order to examine whether 
or not strong naked singularities exist. By the above argument, a light ray that 
emanates from a singularity must lie in the past of the apparent horizon. 

The LTB solution from generic regular initial data results in a shell-focusing 
naked singularity at the center, r = 0. In order to show the existence of a naked 
singularity, we investigate the geodesic equation for an outgoing radial null geodesic 
that emanates from the singularity. For this purpose, we derive the root equation 
that probes the naked singularity, following Joshi and Dwivedi. An outgoing radial 
null geodesic is given as 

(1.34) 

Here we define 

(1-35) 

where a > 1 is determined by requiring x to have a positive finite limit xq as r ^ 0. 
Note that the regular center corresponds to a = 1. Then, from I'Hospital's rule, we 
obtain 

xo = lim — 
r-->o 



1 dR 

lim 



R=xor° 



r— >o ar" ^ dr 

= lim ( R' + e^-^R) . (1-36) 
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Substituting the LTB solution, we obtain 



Xq 



lim 

r— >o ar" ^ 



+ 



R 



VT+f 



(1-37) 



R=xor°' 



In order to obtain the root equation for the LTB solution, we must have an explicit 
expression for R' . By differentiating both sides of Eq. (1-20) with respect to r, we ob- 
tain such an expression of i?', after a straightforward but rather lengthy calculation, 

as 

B! = iV(a;,r)r"-\ 

where N{x, r) is given by 



(1-38) 



N{x,r) = {ri- f3)x + 







7] 



/? ) x^/^G{-Px) 



x 



with 




for fj^O, 
for / = 0, 



l + /3-r? 



(1 +p)l/2^3(a-l)/2 

Therefore, the desired equation is 



+ 



(r?-i/3) G{-p) 

^3(a-l)/2 



a r— >0 



/ + 



VT+f 



(1-39) 

(1-40) 
(1-41) 

(1-42) 
(1-43) 
(1-44) 

(1-45) 
(1-46) 



Note that a should be determined by the requirement that 0{r) have a finite limit 
as r ^ 0. 

For simplicity, we assume that F{r) and /(r) are of the forms 

F(r) = Fgr^ + F^r^ + Fjr"^ + ■■■, (1-47) 
f{r) = f2r^ + hr^ + /er' + • • • • (1-48) 

This implies that the density and specific energy fields arc initially not only finite 
but also analytic at the symmetric center. That is, the initial density and specific 
energy profiles are analytic functions with respect to the locally Cartesian coordi- 
nates. Hereafter we assume F3 > 0, which ensures the positiveness of the central 
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energy density at t = 0. For marginally bound collapse, which is defined by / = 0, 
a positive finite root of Eq. (1-46) is obtained for F5 < as 

with a = 7/3. F5 < means p"{0,0) < 0. Therefore, there exists a naked sin- 
gularity in marginally bound collapse with p"{0.0) < initially. The nakedness of 
the singularity in this spacetime was found numerically by Eardley and Smarr^^^ 
and proved by Christodoulou. For marginally bound collapse with F5 = 0, it is 
easily found that the root equation (1-46) has no positive finite root for any a > 1. 
Therefore, in this case, the singularity is not naked. For homogeneous marginally 
bound collapse, the singularity is not naked, because F5 = 0. For the non-marginally 
bound case / / 0, a similar but more complicated condition for the appearance of 
a naked singularity is obtained, and it was shown that the appearance of a naked 
singularity is generic in the space of this class of functions. is)"20) 

We can also consider a more general class of F{r) and /(r) of the forms 

F(r) = Fsr^' + F^r^ + F^r^ + Fqt'^ + Fyr^ + • • • , (1-50) 
f{r) = f2r^ + fzr^ + Ur^ + hr'' + /gr^ + • • • . (l-Sl) 

This class of the functions corresponds to initial density and specific energy distri- 
butions that are finite but not analytic in general with respect to locally Cartesian 
coordinates. It was shown that a naked singularity also appears from a generic data 
set in this extended space of functions. ^"^^ For marginally bound collapse with 
-F3 > 0, F4 = = and Fq < -(26 + lb^f3)Fl''^ /2, Eq. (1-46) has a finite positive 
root xq with a = 3, and hence the singularity is naked, where xq is given by the root 
of some quadratic equation. 

We can also consider self-similar dust collapse. Self-similarity requires all di- 
mensionless quantities to be functions of f/t for some comoving coordinates {t,f), 
where f is different from r in general. This requirement implies the functional forms 
F{f) = Xf and /(f) = 0, where A is constant. It is found that the singularity is naked 
for < A < 6(26 — 15\/3) and not naked for larger values of A. In fact, through 
some manipulation, it is found that the self-similar case can be classified into the 
marginally bound case with F3 > 0, F4 = = and Fq < 0, which implies that an 
analytic initial density profile is not allowed for the self-similar case. 

With regard to the global visibility of a naked singularity, we can give a simple 
answer. If the Taylor expansions of F(r) and /(r) around the center are such that a 
naked singularity appears, we can immediately construct not only spacetimes with 
globally naked singularities but also those with locally naked singularities by choosing 
the functions F{r) and /(r). In other words, local visibility is determined only by 
the central expansions of the two arbitrary functions, while global visibility depends 
on their functional forms in the whole range < r < 00. In self-similar dust collapse, 
if a singularity is naked, then it is globally naked. 

In contrast to shell-crossing singularities, shell-focusing naked singularities sat- 
isfy the limiting focusing condition for a = 7/3,'^^^ and even both the limiting fo- 
cusing and strong curvature conditions for a = 3 along the first radial null geodesic 
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from the singularity. ^^-^ Irrespective of the value of a, the shell-focusing singu- 
larity satisfies both conditions for time-like geodesies. For the marginally bound 
case, the rcdshift of the first light is finite for a = 7/3^^-* but infinite for a = 3.^^^ 

The LTB spacctimc with a shell-focusing singularity is inextcndible. ^"^^ Detailed 
analysis shows that a naked singularity which may appear in spherical dust collapse is 
ingoing null. This means that there exists a one-parameter family of outgoing 

radial null geodesies that emanate from the singularity, while there exists only one 
ingoing radial null geodesic that terminates at the singularity. Very recently, it was 
shown that nonradial null geodesies, which have nonzero angular momentum, can 
emanate from a singularity if and only if a radial null geodesic emanates from the 
singularity. The appearance of the naked singularity to a distant observer through 
these nonradial null geodesies has been discussed. 

As we have seen above, the collapse of an inhomogeneous dust ball, which 
is given by the LTB solution, results in a shell-focusing naked singularity from 
generic regular initial data. The collapse of a homogeneous dust ball, which is 
called the Oppenheimer-Snyder solution, results in a covered singularity. Though 
the Oppenheimer-Snyder solution was thought to be a typical example of complete 
gravitational collapse, the absence of a naked singularity in this solution turns out 
to be atypical in general spherically symmetric dust collapse. 

1.3. Spherical collapse of realistic matter fields 

It is clear that a dust fluid is not a good matter model, because the effects of 
pressure would not be negligible in actual singularity formation. Here we briefly 
review several examples of gravitational collapse in the presence of the effects of 
some kind of pressure, in the context of singularity formation. We do not describe 
the details of each example, because this is not the aim of this paper. 

1.3.1. Perfect fluid collapse 

The perfect fluid matter model is one of the most natural ways of introducing 
matter pressure. If the pressure is bounded from above, there can appear a shell- 
crossing naked singularity, which is gravitationally weak. We will see below that 
there can appear a strong curvature naked singularity even with unbounded pressure. 

Ori and Piran " '^'^^ investigated the self-similar, spherically symmetric and 
adiabatic gravitational collapse of a perfect fluid with a barotropic equation of state. 
Employing a self-similarity assumption, the equation of state can be restricted to the 
form P = kp, and the Einstein field equations are reduced to a system of ordinary 
differential equations. For < k < 0.4, Ori and Piran found numerically that there 
is a discrete set of self-similar solutions that allow analytic initial data beyond a 
sonic point. These solutions can be labelled by the number of zeroes in the velocity 
field of world lines of constant circumferential radius relative to the fluid element. 
There exists a pure collapse solution among these self-similar solutions, which they 
call the general relativistic Larson-Penston solution. They showed that a central 
naked singularity forms in this Larson-Penston solution for < k ^ 0.0105. They 
also showed that there are naked-singular solutions with oscillations in the velocity 
field for < A: < 0.4. The results of this work were confirmed and extended to 
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< fc < 9/16. This naked singularity is ingoing null. It was shown that this naked 
singularity satisfies both the limiting focusing condition and the strong curvature 
condition for the first null ray. ^^-^ 

However, it is obvious that all initial data sets from which self-similar spacetimes 
develop occupy zero measure in the space of all spherically symmetric regular initial 
data sets. Taking this fact into consideration, there have been discussions that the 
emergence of the naked singularity may be an artifact of the assumption of self- 
similarity. In order to judge the necessity of self-similarity assumption, Harada^^^ 
numerically simulated the spherically symmetric and adiabatic gravitational collapse 
of a perfect fluid with the same equation of state P = kp without this assumption. 
Since null coordinates were used in these simulations, he could detect naked singular- 
ities, not relying upon the absence of an apparent horizon. The result was that naked 
singularities develop from rather generic initial data sets for < fc ^ 0.01, which is 
consistent with the result of Ori and Piran obtained using the self-similarity assump- 
tion. In fact, through further numerical simulations by Harada and Maeda, it was 
found that generic spherical collapse approaches the Larson-Penston self-similar so- 
lution in the region around the center, at least for < < 0.03. This finding is 
supported by a linear stability analysis of self-similar solutions. In other words, the 
Larson-Penston self-similar solution is an attractor solution in the spherical gravi- 
tational collapse of a perfect fluid with the equation of state P = kp, at least for 
< k < 0.03. This work provided the flrst evidence of the generic nature of the 
appearance of naked singularities in spherically symmetric spacetimes with perfect 
fluids. This work also provided the first nontrivial evidence of the attractive nature 
of a self-similarity solution in gravitational collapse. Although the final fate of the 
generic spherical collapse of a perfect fluid with larger values of k is not known, 
Harada analytically showed that the Larson-Penston self-similar solution is no 
longer stable, because of the so-called kink instability, for k > 0.036. 

Undoubtedly, self-similarity plays an important role in certain circumstances of 
gravitational collapse in this attractive case and also in the critical case described 
below. In particular, the former is direct evidence supporting the self-similarity 
hypothesis that spherically symmetric spacetime might naturally evolve from com- 
plicated initial conditions into a self-similar form, which was originally proposed 
in the context of cosmological evolution. (For a recent review of a spherically 
symmetric self-similar systems with perfect fluids, see Ref. 38).) 

1.3.2. Critical collapse 

Critical behavior in gravitational collapse was discovered by Choptuik. "^^^ He 
investigated the threshold between collapse to a black hole and dispersion to infinity 
in a spherically symmetric self-gravitating system of a massless scalar field, i.e., an 
Einstein-Klein-Gordon system. He found so-called critical behavior, such as a scaling 
law for the formed black hole mass, which is analogous to that in statistical physics. 
He also found that there is a discrete self-similar solution that sits at the threshold 
of black hole formation, which is called a "critical solution." Similar phenomena 
have been observed in the collapse of various kinds of matter fields, for example, 
axisymmetric gravitational waves, ^"^^ radiation fluid and more general perfect flu- 
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ids. A renormalization group approach applied by Koike et al. gives a simple 

physical explanation of the critical phenomena and showed that the critical solution 
is characterized by a single unstable mode. Recently, critical behavior was found to 
possibly exist even in the Newtonian collapse of an isothermal gas.^^-* (For a recent 
review of critical phenomena in gravitational collapse, see Ref. 46).) 

A consequence of the mass scaling law for a black hole is the appearance of a 
"zero-mass black hole." A zero-mass black hole can be regarded as a naked singular- 
ity, because the curvature strength on the black hole horizon is proportional to the 
inverse square of the black hole mass. It is obvious that zero-mass black hole forma- 
tion is unstable, because it is a result of exact fine-tuning. On the other hand, if we 
take the limitation of general relativity into consideration, the future predictability 
of classical theory breaks down for a finite-measure set of parameter values in this 
model. The critical collapse of a scalar field provides a very important example of 
naked singularity because it is the first example of a naked singularity that develops 
from regular initial data in the collapse of elementary fields. 

Several theorems regarding the appearance and instability (non-generic nature) 
of naked singularities in a spherically symmetric Einstein-Klein-Gordon system have 
been proved by Christodoulou. 

1.3.3. Collapse of collisionless particles 

One interesting example of spherical collapse is a spherical self-gravitating sys- 
tem of countcrrotating particles, i.e., an Einstein cluster. The static system was 
considered by Einstein, ^'^^ and the corresponding dynamical system was considered 
by Datta, Bondi, and Evans. This system can be constructed by putting 
infinitely many collisionless particles orbiting around the symmetric center with a 
single radial velocity at any given radius, so that the system is spherically symmet- 
ric. Although each particle has conserved angular momentum, the total angular 
momentum vanishes, due to the spherical symmetry. This system is an example of a 
matter field with vanishing radial stress. " The metric functions for a dynamical 
cluster of counterrotating particles are written in terms of an elliptic integral. 
This system has three arbitrary functions, which determine the initial mass distri- 
bution m{r), the energy distribution /(r), and the angular momentum distribution 
L(r). These three are all conserved in this system. It was shown that there appears 
a naked singularity for some class of these arbitrary functions that corresponds to 
regular initial data, although this appearance is not generic in the space of 

all regular initial data sets. In particular, for marginally bound collapse with some 
specific angular momentum distribution m(r) = 4L{r), the metric functions can be 
expressed in terms of elementary functions alone, and this collapse results in naked 
singularity formation, irrespective of the initial density profile. Detailed analysis 
shows that this naked singularity is time-like, unlike naked singularities in spheri- 
cal dust collapse. ^^"^ In fact, this spacetime dynamically asymptotically approaches 
the static model with a central time-like naked singularity. This naked singularity 
satisfies the limiting focusing condition for the first null ray, and even the strong 
curvature condition for a time-like geodesic. 

The above described model is a special realization of a self-gravitating system of 
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collisionless particles, i.e., an Einstein- Vlasov system. This system can be described 
by a distribution function, which obeys the Vlasov equation. (See Ref. 62) for a 
recent review of this system.) The global existence theorem of regular solutions 
with initial data for the distribution function in the Newtonian counterpart of 
this system, the Poisson- Vlasov system, has been proved. This implies that the 
singularity which may form for initial data in the Einstein- Vlasov system is 
not "matter-generated." For a spherically symmetric Einstein- Vlasov system with 
initial data, the global existence theorem of regular solutions with small initial 
data®^) and the regularity theorem away from the symmetric center®^) have been 
proved. 

1.3.4. Other examples 

We should mention several additional examples of naked singularities. Although 
we have restricted our attention to type I matter, there exists a spherical example 
of a naked singularity with type II matter. If we consider imploding "null dust" into 
the center, the spacetime is given by the Vaidya metric. ' For this spacetime, 
it has been shown that naked singularity formation is possible from regular initial 
data. ™) 

There exists a spherical "quasi-spherical" solution in the case of a dust fluid, 
which is called the Szekeres solution. This solution can be regarded as a deforma- 
tion of spherical dust collapse. It has been shown that shell-focusing singularities are 
possible in this solution, and the conditions necessary for the appearance of naked 
singularities are very similar to those in the case of spherical dust collapse. The 
global visibility of this singularity was also examined. 

There have been many analyses with a set of assumptions and the conditions for 
the appearance of a naked singularity are written down in terms of the energy density, 
radial stress and tangential stress. In these analyses, in general the conclusion is that 
naked singularities are possible for generic matter fields that satisfy some energy 
conditions. In this approach, this conclusion is very natural, because there remains 
great freedom in the choice of matter fields, even if some energy condition is imposed. 

Finally, we should mention locally naked singularities in black hole spacetimes. It 
is well known that a Reissner-Nordstrom black hole has a time-like naked singularity 
in the interior of the event horizon. There also exist locally naked singularities in 
the general Kerr-Newmann-de Sitter family of black hole spacetimes. Although these 
singularities do not violate the weaker version of the cosmic censorship conjecture, 
they are inconsistent with the stronger version. For these locally naked singularities, 
there has been a great amount of evidence suggesting that the Cauchy horizon is 
unstable and that it might be replaced by covered null singularities in the presence 
of perturbations. ' 

1.4. Hoop conjecture: axisymmetric and cylindrical collapse 

An important conjecture concerning black hole formation following gravitational 
collapse was made by Thorne. This so-called hoop conjecture states that black 
holes with horizons form when and only when a mass M gets compacted into a region 
whose circumference in every direction in C ^ AnM. He analyzed the causal struc- 
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ture of cylindrically symmetric spacetimes and found remarkably different nature 
from spherically symmetric spacetimes. The event horizon cannot exist in cylindri- 
cally symmetric spacctimc. Therefore, if an infinitely long cylindrical object gravi- 
tationally collapses to a singularity, it becomes a naked singularity, not a black hole 
covered by an event horizon. Then, it is natural to ask what happens if a very long 
but finite object collapses. The hoop conjecture was derived from such a thought 
experiment. It should be noted that the hoop conjecture itself does not assert that 
a naked singularity will not appear. It is a conjecture regarding the necessary and 
sufficient condition for black hole formation. 

As we have reviewed above, many examples of naked singularity formation are 
obtained from the analysis of spherically symmetric spacetime. These examples are 
consistent with the hoop conjecture. The circumference C corresponding to a radius 
R centered at the symmetric center is 2itR. For a central naked singularity of a 
spherically symmetric spacetime, the ratio of C to the gravitational mass M within 
R can be estimated as 



where M^. is a constant and 1 < a < 3. In this case the condition C > inAI coincides 
with the condition that the center is not trapped. Therefore the appearance of a 
central naked singularity in spherical gravitational collapse does not constitute a 
counterexample to the hoop conjecture. 

As for the nonspherically symmetric case, several studies support the hoop con- 
jecture. A pioneering investigation is the numerical simulation of the general rela- 
tivistic collapse of axially symmetric stars. Prom initial conditions that cor- 

respond to a very elongated prolate fluid with sufficiently low thermal energy, the 
fluid collapses with this elongated form maintained. In the numerical simulation, 
evidence of black hole formation was not found. This result suggests that such an 
elongated object might collapse to a naked singularity, and in this case it would 
provide evidence supporting the hoop conjecture. 

Shapiro and Teukolsky numerically studied the evolution of a collisionless gas 
spheroid with fully general relativistic simulations. In their calculations, space- 

times describing collapsing gas spheroid were foliated by using the maximal time 
slicing. They found some evidence that a prolate spheroid with a sufficiently elon- 
gated initial configuration and even with a small angular momentum, might form 
naked singularities. More precisely, they found that when a spheroid is highly pro- 
late, a spindle singularity forms at the pole, where the numerical evolution cannot 
be continued. They also found that the singular region extends outside the matter 
region by showing that the Riemann invariant grows there. Then, they searched for 
trapped surfaces, but found that they do not exist. They considered these results as 
indicating that the spindle might be a naked singularity. However, the absence of 
a trapped surface on their maximal time slicing does not necessarily mean that the 
singularity is indeed naked. To be established whether it is naked or not, we would 
need to investigate the region of the spacetime future of the singularity. 

Wald and lyer^^^ proved that even in Schwarzschild spacetime, it is possible to 
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choose a time slice that comes arbitrarily close to the singularity, yet for which no 
trapped surfaces exists in its past. A simple analytical counterpart of the model 
of the prolate collapse studied numerically by Shapiro and Tcukolsky is provided 
by the Gibbons- Penrose construction. '^^^ This construction considers a thin shell 
of null dust collapsing inward from a past null infinity. Pelath, Tod, and Wald®^) 
gave an explicit example in which trapped surfaces are present on the shell, but 
none exist prior to the last flat slice, thereby explicitly showing that the absence 
of trapped surfaces on a particular, natural slicing does not imply the absence of 
trapped surfaces in spacetime. 

It should be noted that there are problems in the hoop conjecture to prove it as a 
mathematically unambiguous theorem, such as how we should define the mass of the 
object and the length of the hoop. Although we have not found a counterexample to 
the hoop conjecture, it is necessary to obtain a complete solution to these problems. 

§2. Naked singularity formation beyond spherical symmetry 

Many of the examples of naked singularity formation have been obtained by 
the analyses under the assumption of spherical symmetry. Investigation of the non- 
spherical systems have great significance in regard to gravitational wave radiation 
emitted from a forming naked singularity, the stability of the 'nakedness' of a naked 
singularity, and the hoop conjecture. In this section we describe the eflforts made to 
go beyond the case of spherical symmetry. 

2.1. Linear perturbation analyses of the Lemaitre-Tolman-Bondi spacetime 

As mentioned in the previous section, central shell-focusing naked singularities 
should appear from 'generic' initial data in the LTB spacetime. However, there are 
some unrealistic assumptions in this model, e.g., pressureless dust matter, spherical 
symmetry, and so on. Here we consider whether spherical symmetry is essential 
to the appearance of a shell-focusing naked singularity in the LTB spacetime. At 
the same time, we investigate whether a naked singularity can be a strong source 
of gravitational wave bursts. For this purpose, we introduce asphericity into the 
LTB spacetime using a linear perturbation method. The angular dependence of 
perturbations is decomposed into scries of tensorial spherical harmonics. Spherical 
harmonics are called even parity if they have parity ( — 1)' under spatial inversion and 
odd parity if they have parity (—1)^+^. Even and odd perturbations are decoupled 
from each other in the linear perturbation analysis. 

We consider both the odd- and even-parity modes of these perturbations in 
the marginally bound LTB spacetime and examine the stability of the nakedness of 
that naked singularity with respect to those linear perturbations. We also attempt to 
investigate whether the naked singularity is a strong source of gravitational radiation. 

89), 90) 

2.1.1. Odd-parity perturbations 

Here we consider a marginally bound collapse background for simplicity. (See 
§1.2.1 and 1.2.2 for the background solution.) We use gauge invariant perturbation 
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formalism for a general spherically symmetric spacetime established by Gerlach and 
Sengupta. ^-"^^'^^^ Their formalism is described in Appendix A. Hereafter we follow 
the notation and conventions given there. The linearized Einstein equations (A-21) 
and (A-22) can be rewritten as 



■• 1 1 /R' R"\ ( R R'\ ■ 

- = ^ ( 6- - — ) V. - ( 6- + - ) V'. - 4 



i?'2 \ R Rl 



167r _ / r'^p(r)U(r) 



R R' 



R!_R 1 (r: 

RfR^2 [r 



R'R^ 



R'R^ 



(2-1) 



where p{r) = /9(0, r) is the background density profile at t 
the perturbation of the four-velocity as 



5u^ = {d,Q,U{r)S, 



0, and U (r) characterizes 



(2-2) 



where Sa is the odd-parity vector harmonic function. Here the dot and prime denote 
partial derivatives with respect to t and r, respectively. We introduce the gauge 
invariant variable 



i?2 



(2-3) 



where A = e'^ = i?' in the line element (1-1). We solve this partially differential 
equation numerically. 

Let us consider the regularity conditions for the background metric functions 
and gauge-invariant perturbations at r = 0. Hereafter, we restrict ourselves to 
the axisymmetric case, i.e., m = 0. Note that this restriction does not lessen the 
generality of our analysis. Further, we consider only the case in which the spacetime 
is regular before the appearance of the singularity. This means that, before naked 
singularity formation, the metric functions R{t, r) and A(t, r) behave near the center 
according to 



R 
A 



Rc{t)r + 0{r^) 
Rcit) + Oir^). 



(2-4) 
(2-5) 



To investigate the regularity conditions of the gauge- invariant variables ka and Lq, 
we follow Bardeen and Piran. The results are 



ko 
ki 



Le(i)r'+i + 0(r'+3), 



(2-6) 
(2-7) 
(2-8) 



Prom Eqs. (2-3), (2-4), (2-5), (2-7) and (2-8), we find that i/js behaves near the center 
as 



i^s — i^sc{t)r + 0{r^) 



for / > 2, 
for I = 1. 



(2-9) 
(2-10) 
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In the case I > 2, the coefficient ipsc{t) is related to Rc{t) and koc{t) as 

From the above equations, we note that only the quadrupole mode, Z = 2, of ips does 
not vanish at the center. 

We now comment on the behavior of the matter perturbation variable Lq around 
the naked singularity on the slice t = tQ. The regularity conditions for Lq and p 
determine the behavior of U{r) near the center as 

U{r) oc r'+\ (2-12) 

This property does not change even if a central singularity appears. However, the 
r dependence of R and A near the center changes at that time. Assuming a mass 
function F(r) of the form 

F{r) = Fs + F„+3r-"+^ + • • • , (2-13) 

we obtain the relation 

i.W =*o + inr" + --- (2-14) 

from Eqs. (1-16) and (1-29), where n is a positive even integer and to ^ind t„ are 
constants. After substituting this relation into Eq. (1-28), we obtain the behavior of 
R and A around the central singularity as 

R{to,r) ocr^+'i" (2-15) 

and 

^(to,r)ocri" (2-16) 

on the slice t = tQ. As a result, we obtain the r dependence of Lq around the center 
when the naked singularity appears as 

Lo(to,r)ocr'-2"+i. (2-17) 

For example, if / = 2 and n = 2, then Lq is inversely proportional to r and diverges 
toward the central naked singularity. Therefore, the source term of the wave equation 
is expected to have a large magnitude around the naked singularity. Thus the metric 
perturbation variable ipg as well as the matter variable Lq may diverge toward the 

naked singularity. 

In place of the (t, r) coordinate system, we introduce a single-null coordinate 
system, (u, f), where u is an outgoing null coordinate chosen so that it agrees with 
t at the symmetric center, and we choose r = r. We perform the numerical integra- 
tion along two characteristic directions. Therefore we use a double null grid in the 
numerical calculation. In this new coordinate system, {u,f), Eq. (2T) is expressed 
in the form 

^ = _ « U' + ^_RRR' - ^-R'R'] ^,-^\^-1(1-R) 

du R 2 4 ^ 2 R'^ R\ J 
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Sttq f r^p{r)U{r) \ 

df^s = ^(t>s - 3^ (l + R) ^s, (2-19) 

where the ordinary derivative on the left-hand side of Eq. (2-18) and the partial 
derivative on the left-hand side of Eq. (2-19) are given by 

= ^9,-^dr, (2-20) 

9f = -§4^^* + 9r = R'dt + dr. (2-21) 
{dtu)r 

Also, (l)s is defined by Eq. (2-19) and a is given by 

« = (2-22) 

{dtU)r 

We assume ips vanishes on the initial null hypersurface. Therefore, there exist 
initial ingoing waves that offset the waves produced by the source term on the initial 
null hypersurface. In Ref. 89), it is confirmed that this type of initial ingoing waves 
propagate through the dust cloud without net amplification, even when they pass 
through the cloud just before the appearance of the naked singularity. Therefore 
these initial ingoing waves are not significant to the results of this perturbation 
analysis. 

We adopt the initial rest mass density profile 

l + exp(-l^) 

where po, ri and r2 are positive constants and n is a positive even integer. With 
this form, the dust fluid spreads all over the space. However, if r 3> ri,r2, then 
p{r) decreases exponentially, so that the dust cloud is divided into a core part and 
an envelope, which can be considered as essentially the vacuum region. We define a 
core radius as 

rcore = ri + y. (2-24) 

If we set n = 2, there appears a central naked singularity. This singularity becomes 
locally or globally naked, depending on the parameters (/90)^1j^2)- However, if the 
integer n is greater than 2, the final state of the dust cloud is a black hole for any 
values of the parameters. Then, we consider three different density profiles that 
correspond to three types of final states of the dust cloud, globally and locally naked 
singularities and a black hole. The outgoing null coordinate u is chosen so that it 
agrees with the proper time at the symmetric center. Therefore, even if a black 
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hole background is considered, we can analyze the inside of the event horizon. The 
corresponding parameter values are given in Table I. Using this density profile, 
we numerically calculate the total gravitational mass of the dust cloud M. In our 
calculation we adopt the total mass M as the unit of the variables. 



Table I. Parameter values for initial density profiles, power law indices and damped oscillation 

frequencies. 





final state 


Po 


ri 


7-2 


n 


power index 


damped oscillation frequency 


(a) 


globally naked 


f X 10"^ 


0.25 


0.5 


2 


5/3 




(b) 


locally naked 


1 X 10"^ 


0.25 


0.5 


2 


5/3 


0.37+0.089i 


(c) 


black hole 


2 X 10"^ 


2 


0.4 


4 




0.37+0.089i 



The source term of Eq. (2-18), 

^('-)-^(^)' (-) 

is determined by U(r). As mentioned above, the constraints on the functional form 
of U{r) are given by the regularity condition of Lq. Prom Eq. (2-12), U{r) should 
be proportional to r'"*"^ toward the center. We localize the matter perturbation near 
the center to diminish the effects of the initial ingoing waves. Therefore we define 
U{r) such that 



r^p{r)U (r) 




(2-26) 



where Uq and are arbitrary constants. In our numerical calculation we chose r^ to 
be ^core/2- This choice of has no special meaning, and the results of our numerical 
calculations are not sensitive to it. 

First we observe the behavior of tpg at the center. The results are plotted in 
Fig. 1. The initial oscillations correspond to the initial ingoing waves. After these 
oscillations, ipg grows proportionally to (to — t)~^ for the naked singularity cases near 
the formation epoch of the naked singularity. For the case of black hole formation, 
ips exhibits power-law growth in the early part. Later, its slope gradually changes, 
but it grows faster than in the case of a naked singularity. For naked singularity 
cases, the power-law indices 6 are determined by dlji^s/dln(to — t) = (to — t)ips/^s 
locally. The results are shown in Fig. 2. From this figure, we see that the final 
indices are 5/3 for both naked cases. Therefore, the metric perturbations diverge at 
the central naked singularity. 

We also observe the wave form of ipg along the line of constant circumferential 
radius outside the dust cloud. The results are shown in Figs. 3-5. Figure 3 displays 
the wave form of the globally naked case (a). Fig. 4 displays the wave form of the 
locally naked case (b), and Fig. 5 displays the wave form of the black hole case 
(c). The initial oscillations correspond to the initial ingoing waves. In the cases of 
locally naked singularity and black hole formation, damped oscillations dominate the 
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Fig. 1. Plots of tps at the center as a function of the time coordinate t. The soUd curve represents 
the globally naked case (a), the dashed curve represents the locally naked case (b), and the 
dotted curve represents the black hole CEise (c). 




Fig. 2. Plots of the local power indices dln^s /d\n{to — t) = {to — t)tps/^s- The solid curve corre- 
sponds to the globally naked case (a), and the dashed curve corresponds to the locally naked 
case (b). The two curves approach the same value near 5/3. 



gravitational waves. We read the frequencies and damping rates of these damped 
oscillations from Figs. 4 and 5, and express them as the complex frequency 0.37 + 
0.089i for the locally naked singularity and black hole cases. These agree well with 
the fundamental quasi-normal frequency of the quadrupole mode {2Mu! = 0.74734 + 
0.17792i) of a Schwarzschild black hole given by Chandrasekhar and Detweiler. 
In the globally naked singularity case (a), we did not see this damped oscillation 
because of the existence of the Cauchy horizon. In all cases, the gravitational waves 
generated by matter perturbations are at most quasi-normal modes of a black hole 
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Fig. 3. Plots of tpa for the globally naked case (a) at J? = 100. In (a), the first oscillation originates 
from the initial ingoing wave. In (b), wc magnify the the right-hand edge, which is just before 
the Cauchy horizon. The dotted lines represent the time at which the observer at i? = 100 
intersects the Cauchy horizon, which is determined by numerical integration of the null geodesic 
equation from the naked singularity. 



P 





Fig. 4. Plots of tps for the locally naked case (b) aX R= 100. In (a), the first oscillation originates 
from the initial ingoing wave. After this oscillation, the damped oscillation dominates, and this 
part of the wave form is magnified in (b). 



that is generated outside the dust cloud. Therefore intense odd-parity gravitational 
waves would not be produced by inhomogeneous dust cloud collapse. It is thus not 
expected that the central extremely high density region can be observed using this 
mode of gravitational waves. 

We can calculate the radiated power of the gravitational waves and thereby gain 
an understanding of the physical meaning of the gauge-invariant quantities. ' 
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Fig. 5. Plots oiips for the black hole case (c) at R = 100. In (a), the first oscillation originates from 
the initial ingoing wave. After this oscillation, the damped oscillation dominates, as depicted in 
(b). 




Fig. 6. Plots of the radiated power P for the globally naked case (a) at i? = 100. The horizontal 
axis is the outgoing null coordinate u. At the Cauchy horizon, this coordinate has the value mo- 



The radiated power P of the quadrupole mode is given by 

3 



P = 



327r 



-R 



,/2 



(2-27) 



Figure 6 displays the time evolution of the radiated power P. The radiated power 
also has a finite value at the Cauchy horizon. The total energy radiated by odd-parity 
quadrupole gravitational waves during the dust collapse should not diverge. 
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2.1.2. Even-parity perturbations 

The behavior of the even-parity perturbations of the LTB spacetime is investi- 
gated in Ref. 98). There are four gauge-invariant metric variables, kab and k, and 
seven matter variables, Tab, Ta, T^, and T^, where a and b run over and 1. The 
energy density p is perturbed by adding the scalar term 5pY, while the four- velocity 
is perturbed by adding the term 

Su^ = {Vo{x'')Y, Vi{x'')Y, V2{x'^)Ya), (2-28) 

where y is a scalar harmonic function. The normalization for the four-velocity yields 
the relation u'^Sun = 0. This relation implies that Vq vanishes. Then, there are only 
three matter perturbation variables, 

Too = Sp{t,r), (2-29) 

Toi = pVi{t,r), (2-30) 

To = pV2it,r), (2-31) 

as the others vanish: 

Til = Ti = = = 0. (2-32) 

Now we can write down the perturbed Einstein field equations for the background 
LTB spacetime. The resultant linearized Einstein equations are given in Appendix 
C. 

We have derived seven diflFerential equations, (C-l)-(C-7), for seven variables 
(four metric and three matter). The right-hand sides of four of these equations 
vanish. We can obtain the behavior of the metric variables through the integration 
of these equations alone. We transform these equations into more convenient forms. 
From Eq. (A-26), we have 

koo = -^kii. (2-33) 

Using this relation and the equations whose right-hand sides vanish, we obtain evo- 
lution equations for the gauge-invariant metric variables as 

R"R 2R'R' 2R"R 





+ + 




R R R' 

koi, 



2 q' R. R- 2 ( R' R\ , 

koi = -§koi-q', (2-36) 
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where q = k — fcoo- If we solve these three equations with some initial data under 
appropriate boundary conditions, we can follow the full evolution of the metric per- 
turbations. When we substitute these metric perturbations into Eqs. (C-1), (C-2) 
and (C-4), the matter perturbation variables Sp, Vi and V2, respectively, are obtained. 

We can also investigate the evolution of the matter perturbations from the lin- 
earized conservation equations 5{T^^^.^) = 0. They reduce to 



5p\ 1 R^p 



Y ikoi + VM -±V2-k-^{k-q), (2-37) 



p J pR^R' \ R' ' " 'J i?2 

Vi = -l{k'-q'), (2-38) 

V2 = -l{k-q). (2-39) 

Integration of these equations gives us the time evolution of the matter perturbations. 
We can check the consistency of the numerical calculation by comparison of these 
variables and those obtained from Eqs. (C-1), (C-2) and (C-4). 

To constrain the boundary conditions in our numerical calculation, we should 
consider the regularity conditions at the center. These conditions are obtained by 
requiring that all tensor quantities be expandable in nonnegative integer powers of 
locally Cartesian coordinates near the center. We simply quote the results. The 
regularity conditions for the metric perturbations are 

k^koity, q^qoity, koi^koity. (2-40) 

For the matter perturbations, the regularity conditions at the center are 

Sp^Spoity, Vi^Vio{t)r, V2^V2o{ty. (2-41) 

Therefore, all the variables we need to calculate vanish at the center. 

First, we observe the behavior of the metric variables q, k and koi and the Weyl 
scalar, which corresponds to outgoing waves 

^4 = C^^pan>'fh''nPm'' (2-42) 



where 



],-^,0,0), (2.44) 

outside the dust cloud. The results arc plotted in Fig. 7. We can see that the 
metric variables q and feoi and the Weyl scalar ^4 diverge when they approach 
the Cauchy horizon. The asymptotic power indices of these quantities are about 
0.88. On the other hand the metric quantity k does not diverge when the Cauchy 
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Fig. 7. Plots of perturbed variables q, k and fcoi and the Weyl scalar ^4 at constant circumferential 
radius R. The results for _R = 1, i? = 10, R = 100, and R = 200 are plotted. The solid curves 

represent the results for R = 1, the dotted curves for R = 10, the dashed curves for R = 100, 
and the long dashed curves for R = 200. u = uo corresponds to the Cauchy horizon. 



horizon is approached. The energy flux is computed by constructing the Landau- 
Lifshitz pseudo-tensor. We can calculate the radiated power of gravitational waves 
from this. The result is given in Appendix B. For the quadrupole mode, the total 
radiated power becomes 

P = (2-46) 
Stt 

The radiated power of the gravitational waves is proportional to the square of k. 
Therefore it seems unlikely that a system of spherical dust collapse with linear pertur- 
bations can be a strong source of gravitational waves within the linear perturbation 
scheme. However, we should also note that the divergence of the linear perturbation 
variables q, koi and 'P4 implies the breakdown of the linear perturbation scheme. 

Second, we observe the perturbations near the center. The results are plotted 
in Figs. 8 and 9. In these figures we plot the perturbations at t — to = —10^^, 
— 10~^, — lO"^, — 10~^, and 0. Before the formation of the naked singularity, the 
perturbations obey the regularity conditions at the center. Each curve in these 
figures displays this dependence if the radial coordinate is sufficiently small. In 
this region, we can also see that all the variables grow according to power laws as 
functions of the time coordinate along the lines constant r. The asymptotic behavior 
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Fig. 8. Plots of perturbed variables q, k and fcoi near the center. The values for At = to — t = 
10-^10-^10-^ 10-^,0 are plotted. The solid curves represent the results for At = 0, the long 
dashed curves for At = 10^^, the dashed curves for At = 10"^, the dotted curves for At = 10"^, 
and the dotted dashed curves for At = 10"''. 



of perturbations near the central naked singularity is summarized as follows: 

q oc Z\^-2•V^ k oc Z\t-^-V, koi oc Zit-^-V, 

^oczlriV, Vio, At-^\ V20, At-^W 
P 

Here At = tQ — t. On the time slice at At = 0, the perturbations behave as 

jt(xr-°•^^ A;oiocrO-92, 



q oc r 



Sp 



oc r 



-1.4 



Vi oc r°■2^ V2 oc r'--". 



1.3 



On this slice q, k and 6p/p diverge when they approach the central singularity. On 
the other hand, fcoi, Vi and V2 go to zero. 

In the cases of locally naked singularity and black hole formation, we expect 
to observe damped oscillation in the asymptotic region outside the dust cloud, as 
in the odd-parity case. The results are plotted in Fig. 10. These figures show that 
damped oscillations are dominant. We read the frequencies and damping rates of 
these damped oscillations from Fig. 10 and express them as the complex frequen- 
cies 0.36 + 0.096i and 0.36 + 0.093z for locally naked and black hole cases, respec- 
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Fig. 9. Plots of perturbed variables dp, V\ , and V2 near the center. The values for At = 
10-^10-^10-^10-^, and are plotted. The solid curves represent the results for At = 0, 
the long dashed curves for At — 10^^, the dashed curves for At = 10~^, the dotted curves for 
At = 10~^, and the dotted dashed curves for At = 10~*. 



tively. These results agree well with the fundamental quasi-normal frequency of the 
quadrupole mode {2Muj = 0.74734 + 0.17792i). ^4) 

2.1.3. Summary of rclativistic perturbations of spherical dust collapse 

We investigated odd-parity perturbations in §2.1.1. We concluded there that 
the Cauchy horizon is not destroyed by gravitational waves, while the properties of 
a shell-focusing naked central singularity may change, for example, the divergence 
of the magnetic part of the Weyl curvature tensor. 

In §2.1.2, we investigated the behavior of the even-parity perturbations in the 
LTB spacetime. In contrast to the results for the odd-parity mode, the numerical 
analysis for the even-parity perturbations shows that the Cauchy horizon should be 
destroyed by even-parity gravitational radiation. The energy flux of this radiation, 
however, is finite for an observer at constant circumferential radius outside of the 
dust cloud. Therefore inhomogeneous aspherical dust collapse appears unlikely as a 
strong source of gravitational wave bursts. 

The difference between odd and even modes seems to originate from the proper- 
ties of matter perturbations. Odd-parity matter perturbations are produced by the 
rotational motion of the dust cloud, and their evolution decouples from the evolution 
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Fig. 10. Plots of perturbed variables q at constant circumferential radius R = 100 in the locally 

naked singularity and black hole cases, q is normalized with respect to its maximum value, and 
the origin of the time variable is adjusted to coincide with the time when q is maximal. 



of metric perturbations. On the other hand, the even-parity matter perturbations 
contain the radial motion of dust fluid and their evolution couples to the metric 
perturbations. These two modes of odd and even parity, couple to each other when 

we consider second-order perturbations. 

2.2. Estimate of the gravitational radiation from a homogeneous spheroid 

Intuitively, at the formation of a singularity, disturbances of spacetime with 
short wavelength will be created. If there is no event horizon, these disturbances may 
propagate as gravitational radiation, so that a naked singularity may be a strong 
source of short wavelength gravitational radiation. Nakamura, Shibata and Nakao 
have suggested that a naked singularity may emit considerable gravitational wave 
radiation. First, they considered Newtonian prolate dust collapse and estimated the 
amount of the gravitational radiation using a quadrupole formula. 

The shape of the homogeneous prolate spheroid is given by 



The quantities a and c obey 







(2-47) 


obey 






d = 


3M , 1 


(2-48) 


c = 


3M , 1 


(2-49) 


a{e) = 


rfl-^-^%n^+^V 
\ 2e 1 - ej 


(2-50) 


7(e) = 


e2 V 2e 1-e/ 


(2-51) 
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The luminosity of the gravitational waves in the frame of the quadrupole formula is 
given by 

^ = ^p-)^/. = ^((«-f-(c-)Y- (2-52) 
The total amount of energy AE is proportional to 

1= f'~\{a-f-{c-ffdt, (2-53) 

where tc is the time at which a becomes zero, and e > 0. It is easily found that 
the integral / diverges as e — > 0. Therefore, under Newtonian gravity with the 
quadrupole formula, an infinite amount of energy is radiated before the formation 
of a spindle- like singularity. Of course, the collapse of a homogeneous spheroid to a 
singularity cannot be properly described with Newtonian gravity, but this result is 
suggestive. 

To extend this result to a relativistic analysis, they modeled spindle-like naked 

singularity formation in gravitational collapse using a sequence of general relativistic, 
momentarily static initial data for the prolate spheroid. It should be noted that their 
conclusion is the subject of debate. 

2.3. Newtonian analysis on the linear perturbation of spherical dust collapse 

The dynamics of perturbations of the LTB spacetime has been re-analyzed in 
the framework of the Newtonian approximation. In order for a singularity of a 

spherically symmetric spacetime to be naked, "the gravitational potential" 2m/ R 
must be smaller than unity in the neighborhood of the singularity, where m is the 
Misner- Sharp mass function and R is the circumferential radius. A central shell- 
focusing naked singularity of the LTB spacetime satisfies this condition, and further, 
the gravitational potential vanishes even at this singularity. The speed of the dust 
fluid is also much smaller than the speed of light both before and at the time of 
the central shell-focusing naked singularity formation. Therefore, the Newtonian ap- 
proximation seems to be applicable, even though the spacetime curvature diverges 
near the singularity. The advantage of the Newtonian approximation scheme is that 
the dynamics of perturbations of the dust fluid and gravitational waves generated by 
the motion of the dust fluid are estimated separately; the evolution of the perturba- 
tions of the dust fluid is obtained using Newtonian dynamics, and the gravitational 
radiation is obtained using the quadrupole formula. Hence, it is possible to make a 
semi-analytic estimate of the gravitational radiation due to the matter perturbation 
of the LTB spacetime if we adopt the Newtonian approximation. This suggests that 
the Newtonian analysis may be a powerful tool in the analysis of some category of 
naked singularities. However, we should stress that the neighborhood of a naked 
singularity is not Newtonian in an ordinary sense, because there is an indefinitely 
strong tidal force. Thus the Newtonian approximation scheme can be used to de- 
scribe the dynamics of the neighborhood of the naked singularity, but the situation 
as a whole is not Newtonian. 

In this subsection, we consider the gravitational collapse of a spherically sym- 
metric dust fluid in the framework of the Newtonian approximation and show that 
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the Newtonian approximation is valid even at the moment of the formation of a 
central shell-focusing naked singularity if the initial conditions are appropriate, as 
in the case of the example in the previous subsection. 

2.3.1. Eulerian coordinates 

In the Newtonian approximation, the maximal time slicing condition and Eule- 
rian coordinates (for example, the minimal distortion gauge condition) are usually 
adopted. The line element is expressed in the form 

dsl = - (1 + 2^n) dT'^ + dR^ + R^dn^, (2-54) 

where #N is the Newtonian gravitational potential, and we have adopted a polar 
coordinate system as the spatial coordinates. The equations for a spherically sym- 
metric dust fluid and Newtonian gravitational potential <Pn are 

dTP+^dRiR^pV) = 0, (2-55) 
drV + VOrV = -Or^n, (2-56) 
-^dRiR'dR^^) = 47rp, (2-57) 

where V is the velocity of the dust fluid element. The assumptions in the Newtonian 
approximation are 

\V\ < 1 and I^nI < 1, (2-58) 

and further 

\dTV\ <^ \dRVl I^t^nI < |5fl#N| and |Otp| < (2-59) 

2.3.2. Lagrangian coordinates 

For the purpose of following the motion of a dust ball, Lagrangian coordinates 
are more suitable than Eulerian coordinates. The transformation matrix between 
Eulerian and Lagrangian coordinate systems is given by 

dT = dT, (2-60) 
dR = RdT + R'dx, (2-61) 

where r and x are regarded as independent variables, the dot represents a partial 
derivative with respect to r, and the prime represents a partial derivative with respect 
to X. Then the line element in the Lagrangian coordinate system is obtained as 

dsl = - (l + 2^N - R^) dr^ + 2RR'dTdx + R'^dx'^ + R^df2^. (2-62) 

The equations for the dust fluid and Newtonian gravitational potential are 

F{x) = Stt r pR'R^dx, (2-63) 
Jo 

V' = R? = f{x) + ^, (2-64) 

< = ^^(^)' (2-65) 
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where f{x) and F{x) are regarded as arbitrary functions. Since the equation for the 
circumferential radius R is the same as that in the LTB spacetime, its solution in 
the case of the marginally bound collapse, /(x) = 0, has the same functional form 
as Eq. (1-28), 



R = {^^y [TRix) - t]K (2-66) 



where tii(x) is an arbitrary function that determines the time of singularity forma- 
tion. 

Here we consider the Newtonian approximation of the example given in §2.1. 
Therefore we choose the time of the singularity formation as 

™3/2 rj- 

TRix) = (2-67) 

so that R is equal to .x at r = 0. For the initial density configuration, we adopt the 
same functional form as Eq. (2-23), 



where xi and X2 are positive constants. The above choice guarantees the regularity 
of all the variables before singularity formation and that a central shell-focusing 

singTilarity is formed at r = 1. 

Imposing the boundary condition ^ for x ^ oo, the solution of Eq. (2-65) 
can be formally expressed as 

^N = ^Ni(r,x)+^N2(r), (2-69) 

where 

$^^{r,x) = £ ^Fdx, (2-70) 

Here it is worthwhile noting that the right hand side of Eq. (2-65) at x = diverges 
at the time of the central shell-focusing naked singularity formation, 

— ' ^1, x-^/2 for X — ^ at T = 1, (2-72) 

97-7- ' 



where 

_ 1 d?TR{x) 
= 2^^ 



(2-73) 

However, since the power index of x is larger than — 1, itself is finite at x = 0, 
even at the time of the central shell-focusing singularity formation, r = 1. 
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In order for the Newtonian approximation to be successful, temporal derivatives 
of all the quantities should be always smaller than their radial derivatives. Here we 
focus on the neighborhood of the central shell-focusing naked singularity only. For 
this purpose, we introduce a new variable w defined by 



w 



= St-^/^x, (2-74) 



where dr = (1 — T)/Tji(^iy Then, we consider the limit r — 1 as tt; is kept constant. 

It should be noted that x also goes to zero in this limit. The mass function F, 
rest-mass density p, and circumferential radius R behave as 

F^^wV/\ (2-75) 

P ^^Rii)(^ + 7t«')"'(l + wY'Sr-^ (2-76) 

R ^ r^J^^^^il + w'flHr'l^ (2-77) 

All these variables are proportional to powers of (5r, and the coefficients are functions 
of u). It is easy to see that the derivatives of these quantities with respect to r or 
X also have the same basic functional structure with respect to 5r and w. Thus the 
X dependent part, #nij of the Newtonian gravitational potential also behaves in the 
manner 

'Z'Ni — xZ-NiH^rS (2-78) 

where i is a constant and 0ni is a function of w. Substituting the above equation into 
Eq. (2-65) and using the asymptotic behavior described by Eqs.(2-75) and (2-77), we 
obtain 

= ^^^f + ^^') (2-79) 
27rJ/J)(l + w;2)5/3 

In order for the dependences on 5t of two sides of the above equation to agree, i 
must be equal to 1/3. Then, integration of Eq. (2-79) leads to 

^NiH = 0/3,/, ,,,,3 {7(1 + - 9(1 + ^^f"' + 2} . (2-80) 

In order to see the asymptotic dependence of <?n2('7") on r, we differentiate 
Eq. (2-71), obtaining 

The integrand on the right hand side of the above equation behaves near the origin 
at the time of the central shell-focusing singularity formation as 




for X — >0 at r = l. (2-82) 

Therefore, the integral in Eq. (2-81) does not have a finite value at r = 1. Since, 
as shown above, this divergence comes from the irregularity of the integrand at 
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the origin, x = 0, we estimate the contribution to the integral near the origin in 
Eq. (2-81). We again consider the Hmit St ^ diS w is kept constant and obtain 

/ -^\ -^dx — > ^7Tr<^^" ' / 7 ^TTTT = ^7rf>^^^ • (2-83) 

Jo m^lR 9rjg) io (l + «;2)5/3 3^5/3^ ^ ^ 

Substituting the above equation into Eq. (2-81) and integrating it with respect to r, 
we obtain 

^N2 + ^N2(l). (2-84) 

^iJ(i) 

Therefore, in the limit 5t — > 0, with w constant, <Pn can be expressed as 

We now know that in the limit 5t — 0, with w constant, all the variables behave 

as 

Z{t^x) — > z{w)5t^ + constant, (2-86) 

where z{w) is some function of w, and j is a constant. The derivatives of Z with 
respect to T and R can be expressed in terms of derivatives with respect to r and x 
as 

{dTZ)R = {dTT)RZ + {dTx)RZ\ (2-87) 
{dRZ)T = {dRT)TZ + {dRx)TZ'. (2-88) 

From Eqs. (2-60) and (2-61), we find 

{dTT)R = l, {dRT)T = 0, {dTx)R = -^ and {dRx)T = ^,. (2-89) 



Then, we derive the relation 



{dRZ)T z 

Inserting Eq. (2-86) into the above equation, we obtain 

{dRZ)T 3^1/3^(1 + ^2)1/3 \2 ^ ^ 'UwJ j 

(2-91) 

This equation implies that in the limit St ^ with w constant, the inequality 

\{&rZ)R\ < \idRZ)T\ (2-92) 

holds. Therefore, the validity of the order counting of the Newtonian approximation 
is guaranteed even in the neighborhood of the central naked singularity. 
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Here it is worth noting that in the hmit 5t — > with w fixed, |-^| is much larger 
than \Z'\ as seen from 

^^-^(-w- jz^) Sr-y^ oo. (2-93) 

Z' r^(i) \2 ' dzj ^ ' 

Hence, the vicinity of a central shell-focusing naked singularity is not Newtonian in 
the ordinary sense. 

2.3.3. Basic equations of even mode perturbations 

We now consider nonspherical linear perturbations in the system of a spherically 
symmetric dust ball. First, we consider perturbations in the Eulerian coordinate 
system. Here, the line element is written as 

dsl = - (1 + 2#N + 2(5^n) dT^ + dR^ + R'^df2^, (2-94) 

where (5#n is a perturbation of the Newtonian gravitational potential. Using the 
transformation matrix given by Eqs.(2-60) and (2-61), we obtain the perturbed line 
element in the background Lagrangian coordinate system as 

dsl = - (l + 2#N + 25#N - R^) dr^ + 2RR'dTdx + R'^dx'^ + R^df2^. (2-95) 

Hereafter we study the behavior of perturbations in this coordinate system. 
The density p and four-velocity u^^ are written in the forms 

p = p{l + Sp), (2-96) 
u^" = ^1^" + 6u^'. (2-97) 

By definition of the Lagrangian coordinate system, the components of the back- 
ground four-velocity are given by 

(u'^) = (u°,0,0,0) . (2-98) 

Prom the normalization of the four-velocity, we find 

Sv!^ = -(5#N + RR'Su^ . (2-99) 

The order-counting with respect to the expansion parameter e of the Newtonian 
approximation is given by 

6u^ = 0{e^), 6u^ = 0{e), 5p = 0{e^) and 6^N = 0{e^). (2-100) 

Then, the equations for the perturbations are given by 

1 

drSue + de6^N = 0, (2-102) 
■^de (V^^'^^dmS^N) - ^T^pSp = 0, (2-103) 



drSp + -^de [pV^du^] = 0, (2-101) 
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where 

= R'R^smO, (2-104) 

and 7^"* are the contravariant components of the background three-metric. 

Here we focus on axisymmetric even mode of perturbations. Hence the pertur- 
bations we consider are expressed in the forms 

Sp = J2\ii)ir,x)Pi{cose), (2-105) 

I 

S^N = ^A^(l){T,x)Pi{cose), (2-106) 

5ui = J2U,^i){T,x)Piicose), (2-107) 
I 

Su2 = Y,U0^i){T,x)^Pi{cose), (2-108) 

Su3 = 0. (2-109) 
From Eqs. (2-101), (2-102) and (2-103) we obtain 

^P(0 + ^ {^^^m) - + = 0' (2-110) 

{/,(;)+ Z^^(,) = 0, (2-111) 

ile^i) + A^^^i) = 0, (2-112) 

ih^ ^^(o) - + - 4^^"^.(0 = 0- (2-113) 



Comparing the basic equations for the relativistic perturbations with Z = 2 in §2.1.2 
to the above equations, we find the following correspondence between the Newtonian 
and relativistic variables: Z\#(2) = — 2K, Ux(2) = —Vi and Ug(2) = —V2. 

2.3.4. Mass-quadrupole formula 

Hereafter we focus on the quadrupole mode, I = 2, and therefore omit the 
subscript (/) specifying the multipole component of the perturbation variables. The 
mass-quadrupole moment Q^m is given by 

Qem = Jp (xeXm - ^R^6em^ Sx = ^Q(r)diag(-1, -1, 2), (2-114) 

where 

Q{T)= / pApR'^dR. (2-115) 

Jo 

For a function g{T, x) with sufficiently rapid decay as x ^ 00, we find that 

— J^ 9dR = J^ {9rg)RdR = J^ ig-^g'jR'dx (2-116) 

= J (^R'g + R'gjdx-[Rg\^ = J dr{R'g)dx. (2-117) 
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Using the above formula, we obtain 

rinT-n /-oo am i poo f)m 

^ (pApR'R^) dx = — F'^ (a.R') dx. (2-118) 

The power Lqw carried by the gravitational radiation at the future null infinity 
T + — oo is given by 

where u = T — R \s the retarded time. The Wcyl scalar carried by outgoing 
gravitational waves at the future null infinity is estimated as 

= -Cabcdn'^fn'n^fn'' = -^^^ sin^ 9, (2-120) 

where Cabcd is the Weyl tensor, and n" and are two of the null tetrad basis 
vectors, whose components in spherical polar coordinates are given by 

M = ^(1,-1,0,0), (2-121) 

(m^) = ^(0,0, sin 6*). (2-122) 

v2 

It should be noted that the power Lqw is proportional to the square of the third- 
order derivative of Q{u), while the Weyl scalar 1^4 is proportional to the fourth-order 

derivatives of Q{u). 

2.3.5. Asymptotic analysis of the perturbations 

We assume that all of the perturbation variables are regular before the central 
shell-focusing singularity formation and hence can be written as 

Ap = x'^Ap*, (2-123) 

= xU*, (2-124) 

Ue = x'^U;. (2-125) 

A^ = x'^Al, (2-126) 

where each variable with an asterisk is given in the form of a Taylor series with 
respect to z^. 

In order to obtain information concerning the asymptotic behavior of the mass- 
quadrupole moment, we should carefully examine the asymptotic behavior of the 
perturbation variables near the origin. For this purpose, we introduce w defined in 
Eq. (2-74) and then consider the limit r ^ 1 with fixed w. Since all the background 
variables appearing in the equations of the perturbations are proportional to some 
powers of 5t and their coefficients are functions of as given in Eqs. (2-75)-(2-77) 
and (2-85), we expect that the perturbation variables also behave in the same manner 
as the background variables and hence we assume 

Al = S*Jw)St-p, (2-127) 
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Ul = ^d^ix^U;) = r'J^l^ (w^ + 2u;) 5r-^ (2-128) 

U*e=T](^,)U*e{w)5r-\ (2-129) 

(2-130) 
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Now, by virtue of our knowledge about the asymptotic forms (2-127)-(2-130), 
a rigorous analysis about the evolution of the mass-quadrupole moment is possible. 
Substituting Eqs. (2-127)-(2-130) into Eqs. (2-110)-(2-113), and using the asymptotic 
behavior of the background variables (2-75)-(2-77), we obtain 



d5* \ 



18 d \ w'^{l + w'^f'^ d 2 
'^1^ \ (3 + 7u;2)2 d^^"^ ''^> ^ 

I2ul 



U;2(l+u;2)4/3 



rill* \ 

+ 2qu*g St-"-^ + 245r-" = 0, 
dw I 



(2-131) 
(2-132) 



and 



d [i(;2(l+i(;2)5/3 d ^ ^ 



dw 



3 + 7w'^ dw 



{w'd%) 



2i(;2(3 + 7w;2) 

3(l+y;2)l/3 



(2-133) 



Since the powers of St should be balanced in each equation, we obtain 



q=p- 



and 



r = p — 



(2-134) 



Equations (2-131)-(2-133) constitute a closed system of ordinary differential equa- 
tions. 

Through an appropriate manipulation, we obtain a single decoupled equation 
for Uq as 

d'^Un d^Un d^ut du*n ^ ^ , H 

+ + +c^-^ + coul = 0, (2-135) 



dy 



where y = w"^ and 

Co = -2{-24(9 + 26y + 21^^) + 3^^(63 + 4Uy + 1016^^ + 1106?/^ + 44l/) 
+ g(252 + 13232/ + 3149y2 + 3465y3 + 1323y^)} {9y^{l + yf{3 + Tyf}'^ , 

ci = {378 + 21967/ + 4758y^ + 2006y^ - 4424^"^ - 3234y^ 
+ 3g2(l + ?/)^(189 + 1035y + 1911y2 + 1225^^) 
+ g(1323 + 7893y + 18966^^ + 21202^^ + 9247y^ + 441^^)} 
X {18y^{l + yf{S + 7yfy\ 
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C2 = {1269 + 7731y + 18453y^ + 19565/ + 7350/ + 9q^{3 + Wy + 7y^f 
+ 6g(153 + 993?/ + 2387y2 + 2527^^ + 980y^)}{9y\l + yf{2, + lyf}'^, 

C3 = {159 + 506y + 427^^ + 12g(3 + lOy + 7y2)}{6y(l + 2/)(3 + 7^)}"^ 

With appropriate boundary conditions, we can numerically solve Eqs. (2-131)-(2-133) 
as a kind of the eigenvalue problem to obtain q and the solution for Uq. 

The boundary condition at y = for Eq. (2-135) is uniquely determined by the 
Taylor expandability with respect to y in terms of the unknown parameter q and the 
normalization condition UQ\y=Q = 1. We numerically integrate Eq. (2-135) outward 
from y = using the fourth order Runge-Kutta method. The behavior of Ug depends 
on the value of q. 

Considering the behavior of Eq. (2-135) in the limit y ^ oo and imposing the 
condition that Ug is nonzero and finite for < x < e at 5t = 0, where e is a positive 
infinitesimal number, we obtain the outer numerical boundary as 

Ug — > const X y~'^ as y — > oo. (2-136) 

The numerical calculation reveals that the above behavior is realized when 

q = 0.3672 . (2-137) 

From Eq. (2-134), we obtain 

p = 1.701 and r = 1.367 . (2-138) 

The above values agree with the results obtained from the direct numerical simulation 
of partial differential equations (2-110)-(2-113) quite well. 

Now we examine the mass-quadrupole moment Q{T) and its time derivatives 
d"^Q/dT™' . In order to find the contribution of the central singularity to d™'Q/dT'^, 
we consider the integrand on the right hand side of Eq. (2-118). Using Eqs. (2-75), 
(2-77) and (2-127), we obtain 

^From the above equation, we obtain 

I^^Hr,x) = ^ {f'A,R') i_^4/3p^^l3/3-p-n.^26/3-2p-2^ 

|^i;2(p-4/3) ^1 ^ y;2)4/3^*^^^| _ (2-140) 

We consider the integral of J^"*^ from a; = to a; = Xq to determine the con- 
tribution of the central shell-focusing naked singularity to the time derivatives of 



(2-139) 
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the mass-quadrupole moment. Here we take the Umit St with Wg = XqSt 
constant, and then consider the hmit — *■ oo. In this way, we obtain 



jQ Jo 



2 



w 

lo 

2-m 

_^.l/3-m^^29/6-p-m 



3 ^(1) 



The above equation and Eq. (2-138) show that the contribution of the central singu- 
larity to dJ'^Q/dT"^ diverges for r ^ 1 if and only if m is larger than or equal to four. 
This result and the quadrupolc formula imply that the metric perturbation corre- 
sponding to the gravitational radiation and its first-order temporal derivative are 
finite, but the second-order temporal derivative diverges. Hence the power Lqw of 
the gravitational radiation is finite, but the curvature <p4 carried by the gravitational 
waves from the central naked singularity diverges (see Eqs. (2-119) and (2-120)). 
This conclusion agrees with the relativistic perturbation analysis. Further, we find 
that in the limit of r — ^ 1, 

1^^4 = -f ^ oc 6r''^-^ cx (1 - r)-°-«6^ (2-141) 

This result is also consistent with the relativistic perturbation analysis in §2.1.2. 
2.3.6. Summary of Newtonian perturbations of spherical dust collapse 

We analyzed the even-mode perturbations of Z = 2 for spherically symmet- 
ric dust collapse in the framework of the Newtonian approximation and estimated 
the gravitational radiation generated by these perturbations using the quadrupole 
formula. Since we treat separately the dynamics of the matter perturbations and 
the gravitational waves in the wave zone, we can estimate the asymptotic behavior 
semi-analytically, and we obtain the results by solving gentle ordinary differential 
equations. This is the great advantage of the Newtonian approximation. 

^Prom this analysis, we found that the power carried by the gravitational waves 
from the neighborhood of a naked singularity at the symmetric center is finite. How- 
ever, the spacetime curvature associated with the gravitational waves becomes infi- 
nite, in accordance with the power law. This result is consistent with the relativistic 
perturbation analysis in §2.1.2. Furthermore, the power index obtained from the 
Newtonian analysis also agrees with that obtained from the relativistic perturbation 
analysis quite well. 

The agreement between the results of the Newtonian and relativistic analyses 
suggests that the perturbations themselves are always confined within the range to 
which the Newtonian approximation is applicable. Here we focus on the metric per- 
turbation, A*^. Since the asymptotic solution of has the same form as Eq. (2-86), 
we immediately find that in the limit (5r — ^ with w constant. 



(2-142) 
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and hence the assumption ^ of the Newtonian approximation is vahd 

in the Eulerian coordinate system. We can also find the second order derivatives. In 
the same hmit, we find 



d\A% dndrAl 
From the above equations, we obtain 

dlA% 



o.St'/" and <x6t'/^ (2-143) 



r)'^ A* 



OC 



6t^/^. (2-144) 



The above equation implies that in the limit of (5r ^ with fixed w, the inequality 

\d^A%\ « \dlA*^\ (2-145) 

is also satisfied. This inequality implies that the wave equation for the metric pertur- 
bation A^ is approximated well by a Poisson type equation if we adopt the Eulerian 
coordinate system. 

However, the gravitational collapse producing the shell-focusing globally naked 
singularity is not Newtonian in the ordinary sense. The same is true for the perturba- 
tion variables, because \A4,/A'^\ » 1 in the limit St ^ with fixed w. Even though 
the Newtonian approximation is valid in the Eulerian coordinate system, Newtonian 
order counting breaks down if we adopt the Lagrangian coordinate system as the 
spatial coordinates. 

2.4. Cylindrical collapse and gravitational radiation 

It has long been known that collapsing cylindrically symmetric fluids form naked 
singularities. Such examples are not considered as direct counterexamples to the 
cosmic censorship conjecture, because these spacetimes are not asymptotically fiat. 
However, there is an expectation in which the local behavior of prolate collapse to a 
spindle singularity is very similar to that of infinite cylindrical collapse. For this rea- 
son, properties of cylindrical collapse have been studied in this context. Apostolatos 
and Thornc ^'^^^ investigated the collapse of a counterrotating dust shell cylinder and 
showed that rotation, even if it is infinitesimally small, can halt the gravitational 
collapse of the cylinder. Echeveria ^''^^ studied the evolution of a cylindrical dust 
shell analytically at late times and numerically for all times. It was found that the 
shell collapses to form a strong singularity in finite proper time. The numerical re- 
sults showed that a sharp burst of gravitational waves is emitted by the shell just 
before the singularity forms. Chiba showed that the maximal time slicing never 
possesses the singularity avoidance property in cylindrically symmetric spacetimes 
and proposed a new time slicing that may be suitable to investigate the formation 
of cylindrical singularities. He numerically investigated cylindrical dust collapse to 
elucidate the role of gravitational waves and found that a negligible amount of grav- 
itational waves is emitted during the free fall time. There seems to be a discrepancy 
in the results of Echcvcria and Chiba with regard to the emission of gravitational 
waves. The origin of this discrepancy may be the difference between the matter 
fields: a thin, massive dust shell in the former case and a dust cloud in the latter. 
However, further careful analysis is necessary. 
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§3. Quantum particle creation from a forming naked singularity 

As we have seen in §1, naked singularities are formed in some models of gravita- 
tional collapse. If this is the case for more realistic situations, then what happens? 
The existence of naked singularities implies that the high-curvature region due to 
strong gravity is exposed to us. If naked singularities emit anything, we may obtain 
information of some features of quantum gravity. In this context, particle creation 
due to effects of quantum fields in curved space will be one of the interesting possibil- 
ities. Using a semi-classical theory of quantum fields in curved space, Hawking ^^^^ 
derived black body radiation from a black hole formed in complete gravitational col- 
lapse. Ford and Parker ^'^^^ calculated quantum emission from a shell-crossing naked 
singularity and obtained a finite amount of flux. Hiscock, Williams and Eardlcy ^^^^ 
considered a shell-focusing naked singularity which results from a self-similar implo- 
sion of null dust and obtained diverging flux. Here, we will apply a semi-classical 
theory of quantum flelds in curved space to the collapse of a dust ball, i.e., the LTB 
solution. 

3.1. Particle creation by a collapsing body 
3.1.1. Power, energy and spectrum 

We consider both minimally and conformally coupled massless scalar fields in a 
four-dimensional spacetime which is spherically symmetric and asymptotically flat. 
Let r, R, 9, (f) denote the usual quasi-Minkowskian time and spherical coordinates, 
which arc asymptotically related to null coordinates u and v through u ^ T — R 
and V K T + R. If the exterior region is vacuum and spherically symmetric, it is 
described by the Schwarzschild metric, which is given by 



with R^ = R + 2Mln[{R/2M) - 1]. 

An incoming null ray v = const, originating from past null infinity X~, propa- 
gates through the center becoming an outgoing null ray u = const, and arriving on 
future null infinity at a value u = F(v). Conversely, we can trace a null ray from 
u on T"*" to V = G{u) on J^, where G is the inverse of F. Here, we assume that 
the geometrical optics approximation is valid. The geometrical optics approximation 
implies that the trajectories of the null rays give surfaces of constant phase. Then, 
in the asymptotic region, the mode function which contains an ingoing mode of the 
standard form on T~ is the following 







1 




e-'^''^^^)YUe,ct>) 



(3-4) 



■s/AttujR 
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where we have imposed the reflection symmetry condition at the center. In the 
asymptotic region, the mode function which contains an outgoing mode of the stan- 
dard form on X"*" is the following 



-iu>F(v) _ 



)Ylm{0,4>). 



Note that in the above we have normalized the mode functions as 



(3-5) 



(3-6) 



where the inner product is defined by integration on the space-like hypersurface S 
as 

if,9) = -i [ U9l-f,t.9*)dS^. (3-7) 
Using the above mode functions we can express the scalar field (f) as 

^ = E / dj^'WlmK'lm + <L<nm). (3-8) 
l,m 

<P = T.J d^'ia^jLuZ'^L + aTiL<'i:n)- (3-9) 



l.m ' 



According to the usual procedure of canonical quantization, we obtain the following 
commutation relations 



"wimp "■uj'l'm' 

out outt 
"o/Zm' "'ui'l'm' 



5{uj - uj')6ii'6„ 

6{UJ - J)5iii5r, 



(3-10) 
(3-11) 



where it is noted that the Lagrangian in the Minkowski spacetime is common for both 
minimally and conformally coupled scalar fields. Here, aj^^,„j and a°'^*^ arc interpreted 
as annihilation operators corresponding to in and out modes, respectively. Then we 
set the initial quantum state to in vacuum, i.e.. 



«L"m|0) = 0. 



(3-12) 



The radiated power for fixed I and m is given by estimating the expectation value 
of stress-energy tensor through the point-splitting regularization in a fiat spacetime 

agios) 



Im 



R^dQ 



1 

24^ 



3 fG" 



2\G' 



1 



487r V G' 



G" 



for a minimally coupled scalar field, and 



247r \G' ) 

(3-13) 



(3-14) 
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for a conformally coupled scalar field. Here the prime denotes the differentiation 
with respect to the argument of the function. It implies that the amount of the 
power depends on the way of coupling of the scalar field with gravity. However, if 



G' 



b 



(3-15) 



G' 

holds, the radiated energy from u = a to u = boi a minimally coupled field 

Eim= t Pimdu (3-16) 
J a 

and that of a conformally coupled field 

Elm^ f Plradu (3-17) 
J a 

coincide exactly. The actual power is given by summation of all {l,m). The simple 
summation diverges. This is because we have neglected the back scattering effect by 
the curvature potential which will reduce the radiated fiux considerably for larger I. 
Therefore we should recognize that the above expressions for the power, (3-13) and 
(3- 14), are a good approximation only for smaller I. Hereafter we omit the suffixes I 
and m. 

The spectrum of radiation is derived from the Bogoliubov coefficients which 
relate in and out modes given as: ^^"^^ 

1 // j/ roo 

ao^'c = W, ) = 7r\ - dve'^^^-')-'^ ^ (3-18) 

Po^'u. = -(C,«r*) = -7r\ - / t^^e-*"''^'^"*"'"- (3-19) 
The expectation value N{u) of the particle number of a frequency u on J+ is obtained 

by 

roo 

Niu) = / dw'lP^^^f. (3-20) 
Jo 

It is noted that these results are free of ambiguity coming from local curvature 
because the regularization is done only in a fiat spacetime. 

3.1.2. Quantum stress-energy tensor in a two-dimensional spacetime 

We can estimate the vacuum expectation value of stress-energy tensor in a two- 
dimensional spacetime without ambiguity which may come from local curvature in 
contrast to a four-dimensional case. Unfortunately, this is not the case for four- 
dimensional spacetime. For simplicity, wc consider a minimally coupled scalar field 
as a quantum field, although the situation would not be changed for other massless 
fields. 

It is known that any two-dimensional spacetime is conformally flat. Then its 
metric can be expressed by double null coordinates {u, v) as 

ds"^ = -G^{u,v)dudv. (3-21) 
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If the initial quantum state is set to the vacuum state in the Minkowski spacetime 
ds^ = —dudv, then the vacuum expectation value of the stress-energy tensor of the 
scalar field is given by 

(n<.) = ~C{C-')^^, (3.22) 

= — , (3.24) 

where TZ is the two-dimensional scalar curvature. 

As usual, we require that the regular center is given hy u = v and that v coincides 
with the standard Eddington-Finkelstein advanced time coordinate v. We introduce 
the internal double null coordinates U and V. Assuming the relation between internal 
and external null coordinates, 

U = a{u), (3-25) 
V = P{V), (3-26) 

we can obtain all the coordinate relations as 

v = v = P{V), (3-27) 
u = p{U) = p{a{u)). (3-28) 

In terms of the interior double null coordinates, the line element in the interior is 
given by 

ds2 ^ -A^{U, V)dUdV, (3-29) 

and wc further require U = V at the regular center. Then we can transform the 
quantum stress-energy tensor given by {u, v) coordinates to the ones in the interior 
and the exterior double null coordinates. 

Using the coordinate relation equation (3-28) we can transform Eqs. (3-22)- 
(3-24) to the components expressed in the interior coordinates and in the exterior 
coordinates. The components in the interior coordinates are given by 

(Tuu) = FuiP') - Fu{A^) (3-30) 
{Tvv)=FviP')-FviA^) (3-31) 

{Tuv) = -^{lnA''),uv (3-32) 

where 

The P' in Eq. (3-30) should be considered as p'{U). 

In the exterior vacuum region, the line element is given by 

ds'^ = - (l - ^) dudv, (3-34) 
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where u and v are the Eddington-Finkelstein null coordinates. Then, we obtain 

(Tun) = -Fu{D^) + a'^Fu{(3') + F^{a') (3-35) 
(r,„) = -F„{D^) (3-36) 

{Tuv) = -^{InD^),^^, (3-37) 

where D'^ = 1 - 2M/R. The /?' in Eq. (3-35) should be also considered as p'{U). In 
this region we obtain 

(lnl^V = -(2|^-^), (3-38) 

F„(Z^=) = F„(0^) = -— (^-^-^j. (3 39) 

It is found that the outgoing part of flux (T„u) has terms which are dependent on 
the internal structure of the collapsing body. 

3.2. Particle creation in analytic dust collapse 
3.2.1. Mapping 

We can find the function G or F by solving the trajectories of outgoing and 
ingoing null rays in the dust cloud and determining the retarded time u and the 
advanced time v through Eqs. (3-2) and (3-3) at the time when the outgoing and 
ingoing null rays reach the surface boundary, respectively. The LTB solution for 
marginally bound collapse is given by Eq. (1-28). The trajectories of null rays in the 
dust cloud are given by the following ordinary differential equation 

I = ±Rr, (3-40) 

where the upper and lower signs denote outgoing and ingoing null rays, respectively. 

We have solved the ordinary differential equation (3-40) numerically. ^^^^ 
We have Tiscd the Runge-Kutta method of the fourth order. We have executed the 
quadruple precision calculations. 

We have chosen the mass function as 

F{r) = Fsr^ + Fgr^ (3-41) 

We have found that the central singularity has been globally naked for very small 
ri, if we have fixed the value of F3 and F5. Although we have calculated several 
models, we only display the numerical results for the model with F^ = 1, F^ = —2 
and Tft = 0.02 in an arbitrary unit because the features have been the same if there 
has been globally naked singularity in the model. The total gravitational mass M 
is given by M = 3.9968 x 10~^ for this model. See Fig. 11 for trajectories of null 
geodesies. We also indicate the location of the singularity, apparent horizon and 
Cauchy horizon in this figure. 

See Fig. 12. In Fig. 12(a), the relation between u and v, i.e., the function G{u) is 
shown, where uq {vq) is defined as the retarded (advanced) time of the earliest light 
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Fig. 11. Null rays interior of the dust cloud in the analytic LTB spacctirnc. Parameters are chosen 
as Fz = 1, = —2 and Vh = 0.02. See text for definitions of these parameters. 



ray which originates from (terminates at) the singularity. In Fig. 12(b), the first 
derivative G'{u) is shown. This imphes that G'{u) does not diverge but converges to 
some positive constant A with < ^ < 1. In Fig. 12(c), it is found that the second 
derivative G"{u) does diverge as u uq. This figure shows that the behaviors of 
the growth of G"{u) are different from each other during early times (10~^ ^uq — u) 
and during late times (0 < — "u ^ 10~^). During late times, the dependence is 
written as 

G" -{UQ - u)-^/'^ . (3-42) 

It is noted that this singular behavior was re-derived analytically under several as- 
sumptions. ^^^^ 

Here we determine the amplitude of the power by physical considerations. We 
assume that the particle creation during late times is dtic to the growth of the cen- 
tral curvature. Then, we consider a constant of proportion which will appear in 
Eq. (3-42). We should note that the coefficient must be written by initial data 
because it must not depend on time. Equation (3-42) demands a dimensionful con- 
stant of proportion. For late times, the only possible quantity is (^o/^o)' where Iq 
the free-fall time given by 

to = \f^"^ (3-43) 
and Zo denotes the scale of inhomogeneity defined as 




(3-44) 



This is because only (to/^o) does not depend on the origin of time coordinate and 
can be formed by local quantities at the center. See Appendix D for the derivation 
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Fig. 12. Plots of (a) G{u), (b) G'{u) and (c) G"{u) for the analytic LTB spacetime. Parameters 
are the same as those in Fig. 11 



and the physical meaning of this constant. Thus we can determine the coefficient up 
a numerical factor as 



G' ^A, 



-1/2 



(3-45) 
(3-46) 



for late times, where < ^4 < 1 and / is a dimensionless positive constant of order 
unity. It implies that there is characteristic frequency of singularity which is defined 
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Fig. 13. Power for minimally and couformally coupled scalar fields in the analytic LTB spacetime. 
Parameters are the same as those in Fig. 11 



as 

u^s ^ 1 (3-47) 
'^o 

The late-time behavior is good for < tto — « <C The above estimate shows a 
good agreement with numerical results. 

3.2.2. Power and energy 

Since G" and therefore G'" diverge as n — > mq, the power of radiation diverges 
for both minimally and conformally coupled scalar fields. For this model, we have 
numerically calculated the power of radiation by Eqs. (3-13) and (3-14). The results 
are displayed in Fig. 13. 

Based on the above physical considerations, we can obtain the formula for the 
power by the particle creation using Eqs. (3-13) and (3-14). As seen in Fig. 12(b), 
we find 

G' « A, (3-48) 
during late times. The power during late times is obtained as 

P^^M/'iuo-u)-^/', (3-49) 

P^^fu;s{uo-u)-\ (3-50) 

Therefore the power diverges to positive infinity for both minimally and conformally 
coupled scalar fields as u ^ uq. The radiated energy is estimated by integration. 
During late times, the totally radiated energy is estimated as 

E « ^M^H^o - u)-''\ (3-51) 
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4o7r iOs{uQ — nj 

Therefore, the totally radiated energy diverges to positive infinity for both minimally 
and conformally coupled scalar fields as uq- However, in realistic situations, we 

may assume that the naked singularity formation is prevented by some mechanism 
and that the quantum particle creation is ceased at the time At. In 

other words, G"/G' tends to vanish for uq — u < At. In such situations, the total 
derivative term in the expression of the power of a minimally coupled scalar field 
gives no contribution to the totally radiated energy. Therefore, the total energy for 
a minimally coupled scalar field and that for a conformally coupled one coincide 
exactly, i.e., 

E = E -^fusln^—. (3-53) 
487r'' LUsAt ^ ' 

3.2.3. Spectrum 

From the above numerical results and physical considerations, we obtain the 
following formula for the function G{u) or F{v): 

G{u) « A{u - uo) - ^Afiul/^{uo - ufl"^ + vq (3-54) 

or 

F{v) « A~\v - vo) + yul'^[A-\vo - v)f'^ + (3-55) 

for late times. 

Now that we have obtained the function F{v), we can calculate the spectrum 
of radiation by Eqs. (3-18), (3-19) and (3-20). See Harada, Iguchi and Nakao^^'^) for 
details. In order to determine the spectrum numerically, we introduce the Gaussian 
window function in the Fourier transformation as usual. 

Although we could use the calculated data for the function F(v), we have used 
the analytic formula which have been derived based on the numerical results for 
convenience of retaining accuracy. Since the diverging power is associated with the 
latc-timc behavior, wc concentrate on the latc-timc radiation. In order to obtain the 
totally integrated spectrum, we extrapolate the function F linearly as 

A''^{v — Vo) + Uq {vo<v) 



where vi is given as 



vo-vi = 



1 ^2 



(3-56) 



AluJ\ (3-57) 



We should note that Eq. (3-56) can be used for uq — u ^ ^7^- I^i f^ct, it turns out 
that we only have to pay attention to the spectrum above Ug. Wc should note that 
there is no radiation for u > uq in this extrapolation. Because the second derivative 
of F is divergent at v = vq, we can only require that the first derivative of F should 
be continuous a,t v = vq. If we allow discontinuity of the first derivative, radiation 
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(a) Totally integrated spectrum 
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Fig. 14. (a) Totally integrated and (b) momentary spectra for A = 0.8 and / = 1. The latter is for 
the moment uo — u = W~^ojJ^. 



due to this discontinuity dominates the radiated energy flux, which is out of our 
concern. 

The obtained spectrum is shown in Fig. 14(a). The parameters are fixed as 
A = 0.8 and / = 1. We should note that the contribution to the total radiated 
energy mainly comes from u ^ lOs- In this figure, we find 



dN 

LO— — OC LO 

duo 



-1 



Therefore, the total energy, which will be obtained by 

dN 



f°° , dN 
E = I dioio— — , 
Jo dcj 



(3-58) 



(3-59) 



is logarithmically divergent, which is consistent with the divergent radiated energy 
obtained based on the point-splitting regularization. 

It is important in estimating the validity of the geometrical optics approximation 
to examine which frequency band dominates the power at each moment. For this 
purpose, we determine the momentary spectrum by a wavelet analysis. We use the 
Gabor wavelet in place of e*^ because it has a very clear physical meaning as the 
Gaussian wave packet. The Gabor wavelet is written as 

V'(x) = ^^e^^e-^'/'^', (3-60) 



47r(7 

where we have set £7 = 8. The result below is not so sensitive to the value of a. 

We have calculated numerically the momentary spectrum. The result is shown in 
Fig. 14(b). This figure displays the contribution of each logarithmic bin of frequency 
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to the power at the moment. The parameters are fixed as ^ = 0.8 and / = 1. We 
have chosen the time of observation uq — u = 10-^a;7^ In this fi ffure, it is found 
that the contribution to the power is dominated by the frequency uo ~ 27r(uo — u)"^. 

3.2.4. Null geodesies and redshifts 

Here we analytically examine null geodesies and redshifts in the LTB spacetime 
and derive some useful results to comprehend the problem. We assume the function 
F{r) of the form Eq. (1-47) with F3 > and F5 < 0. We define 77 as 

ri = t-to. (3-61) 

Moreover, r]s{r) = ts{r) — to and rjahir) = tah{r) — io are approximated for r/^o ^ 1 
as 

(3-62) 



r]s{r) 


4 


f r 


to 






riahir) 




( r 


to 







3/ \toJ \loJ 
For a while, we assume 

^ < 1. (3-64) 

'0 

This implies that we can safely expand in powers of v/Iq and take the leading order 
term in the whole of the cloud. 

Then we examine trajectories of null geodesies in this spacetime. First we con- 
sider the region where 



to 



(3-65) 



is satisfied. This region is approximately recognized as the Friedmann universe. 
From Eq. (3-62), it is found that this condition can be satisfied only for r] < 0. In 
this case, R^r is approximated as 

R^^(^y^\ (3-66) 



to / 

The ordinary differential equation (3-40) can be easily integrated to 

lfi»VfUc.,. (3.67) 



to J 3 V^o/ \ No- 

where Ca± is a constant of integration. It is easily found that each null ray can be 
drawn by a parallel transport of the Ca± = curve along the r-axis direction. 
Next we consider the region where 



2 ^ 
> — 
to 



(3-68) 



is satisfied. This region is not at all approximated by the Friedmann model. In this 
case, R^r is approximated as 

Rr-'T^oj^ r) ■ (3-69) 



22/33 \loJ 
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Fig. 15. Schematic figure of null geodesies around naked singularity in the analytic LTB spacetime 



where Cb± is a constant of integration. It is easily found that each nuh ray can be 
drawn by a parallel transport of the Cs± = curve along the 77-axis direction. 

See Fig. 15, which illustrates the spacetime around (r/,r) = (0,0). Condition 
(3-65) is satisfied in region A, while condition (3-68) is satisfied in region B. The 
boundary of regions A and B will be described by 



where 7 is a constant of order unity. We denote this boundary curve as C . It may 
be kept on mind that this treatment is rather simple. However we believe that the 
present approximation that we divide the spacetime to two regions will be enough 
to comprehend an essence of the problem. 

Then, we concentrate on the behavior of null geodesies around the naked singu- 
larity. See Fig. 15 for the trajectory of null geodesies. From the regular center r/ < 
and r = 0, the outgoing null ray runs region A which is described by the upper sign 
of Eq. (3-67) with Ca+ > 0. Then this null ray goes through boundary C. After 
that the null ray goes into region B and the trajectory is described by the upper sign 
of Eq. (3-70) with Cb+ > 0. Since all null rays with (7^+ > in region B emanate 
from the regular center, the outgoing null geodesic with Cb^ = in region B is the 
Cauchy horizon, i.e., the earliest outgoing null ray from the singularity at r = 0. 
We denote this outgoing null geodesic as ra_. There are infinitely many outgoing 
null geodesies later than n_ which emanates from the singularity at r = 0. In fact, 
these null geodesies form one parameter family. These null rays cannot asymptote 



Then, Eq. (3-40) is integrated to 




(3-70) 




(3-71) 
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to n_ as r — > because these null rays are obtained by a parallel transport of n_ 
in ry-axis direction in region B. Instead, these null rays asymptote to the location of 
the apparent horizon in an approach to the singularity. 

It is clear that the ingoing null ray with Cb- = terminates at the singularity 
r = 0. We denote this ingoing null ray as n+. For small positive Cb-, the ingoing null 
ray in region B crosses boundary C and the null ray in region A becomes described 
by the lower sign of Eq. (3-67) with Ca- > 0, and then terminates at the regular 
center. For small negative Cb-, the ingoing null ray in region B terminates at the 
space-like singularity r > 0. 

Using the above results, we can derive the condition for globally naked singular- 
ity. If n_ reach then the singularity is globally naked, otherwise the singularity 
is locally naked. Noting that the intersection of the apparent horizon with the cloud 
surface is on the event horizon, we find that the condition for the singularity to be 
globally naked is given as 

2 

> 2^/^. (3-72) 

Since — (16/27)x^ is a decreasing function of x > 0, we can find the condition 
for being globally naked as 

^< 0.66(1)'. (3-73) 

This implies that the singularity is globally naked for sufficiently small r^ when and 
^0 are fixed. Thus, it turns out that the present assumption rh/lo <^ 1 is relevant for 
the globally naked singularity. We can translate the above condition to the condition 
for M and Ug as 

M < 6.4 X lO'^uj-^. (3-74) 

It implies that the mass of the dust cloud with globally naked singularity is bounded 

from above. 

Since it is necessary for later discussions, we turn our attention to estimate of 
redshifts. Let be the tangent vector of a radial null geodesic. The frequency a) 
observed by an observer comoving with a fluid element is calculated from k^^ as 

u = -ki^Uf, = -k\ (3-75) 

In region A, null geodesic equations can be integrated as 

k\-r]f'^ « const, (3-76) 

in the lowest order. Therefore a particle is blueshifted for the comoving observer in 
region A. In region B, null geodesic equations can be integrated as 

k* const, (3-77) 

in the lowest order. Therefore a particle is neither redshifted nor blueshifted in region 
B. 
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In the external Schwarzschild spacetime, in a similar way, we obtain the tangent 
vector of a radial null geodesic as 

(i _ = const, (3-78) 



R J 

A static distant observer observes the frequency 

Lo = -k^{T,oo). (3-79) 

Using the matching condition (1-31), the observed frequency u = —k^u^ by the 
comoving observer at the surface is written using the above u as 

(b = (3-80) 

1 ^ a 

where a is defined by 



a^y^, (3-81) 

and R is the circumferential radius of the dust surface when the null ray crosses the 
surface. 

Let a± be the value of a for n±. Let ri± be the value of rj when n± crosses the 
dust surface. Then, from Eq. (1-31), it is easy to see that the Taylor expansion is 
valid around ij = rj+ for v as a function of ij as 

v = vo + T^—{ri-r]+) + 0{{ri-ri+f). (3-82) 
1 + a+ 

For a- < 1, from Eq. (1-31), it is also easy to see that the Taylor expansion is valid 
around = ry- for n as a function of r/ as 

u = uo + -^{v-V-) + Oi{v-V-?), (3-83) 
1 — a- 

or equivalently, 

r/- - r/ = (1 - a-){uo - u) + O{{uo - uf). (3-84) 

See Fig. 15 again. The light ray which is close to n± enters the dust cloud in 
region B (a), enters region A (b), reaches the center (c), enters again region B (d) and 
leaves the dust cloud (e). In order to determine G or F, we need to relate the time 
coordinate tja at which the incoming light ray enters the dust cloud with rje at which 
the outgoing light ray leaves the dust cloud. Within the present approximation, we 
find 

r/_ - ??e ~ ??+ - ??a, (3-85) 

which can be derived using the solutions of null geodesic equations obtained above. 
Then, we find for a_ < 1 

1 - a_ 

v = G{u)^vo-— {uq-u), (3-86) 

1 + a+ 
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c2;>27r|i?-.-.|V2 (3-89) 



or 

u = F{v)^uo-l^^ivo-v). (3-87) 
1 — a_ 

Prom the above discussions, we can determine A = limu_^.uo G' as 

A = ^—^. (3-88) 

3.2.5. Validity of geometrical optics approximation 

The geometrical optics approximation is exact for a two-dimensional spacetime, 
where the line element is given by setting df2^ = in the four-dimensional line 
element (1-1). Therefore, the obtained power, energy and spectrum are exact in this 
sense for the two-dimensional version of the LTB spacetime. 

In a four-dimensional case, the geometrical optics approximation is only an ap- 
proximation because of the existence of curvature potential. The geometrical optics 
approximation is valid for waves with the wave length shorter than the curvature 
radius of the spacetime geometry. In other words, the geometrical optics approxi- 
mation is good if condition 

is satisfied, where the hat denotes the components in the local inertial tetrad frame. 
In region A, nonvanishing components of the Riemann tensor are approximated 

as 

l^trtfl ~ \^me\ ~ l^rOfel ~ \^e^04,\ ~ (-^)~^- (3-90) 
In region B, they are approximated as 

We have already seen that, in region B, cb is kept constant approximately. In region 
A, Cj may be considerably blueshifted. However, we can see that the null ray which 
goes into the dust cloud enters region A and get out of region A at the same time in 
the lowest order, i.e, rji, ^ rjc ^ rid- This implies that the frequency of such a particle 
is kept almost constant all over the dust cloud. The Riemann tensor along the light 
ray reaches the maximum when the light ray goes through region A. Now we can 
write down the condition for the geometric optics approximation to be valid as 

u > 2Tr{-r]c)-\ (3-92) 

This condition can be rewritten by the quantities on J+ using Eq. (3-80) and (3-83). 
The result is 

CO ^ uJcriu) = 2n{uo — u)~^ (3-93) 

for a_ < 1. 

As we have seen, the power at each moment is mainly carried by particles of 
which frequency is about 27r(no — u)~^. Since the transmission coefficient will 
be a function of {ui/ucr)-, it is natural to estimate /^(a; Ucr) as of order unity. It 
implies that the calculations based on the geometrical optics approximation will be 
valid for a rough order estimate. 
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3.3. Particle creation in self-similar dust collapse 
3.3.1. Null coordinates 

As we have seen in §1.2, the self-similar dust collapse corresponds to marginally 
bound collapse with F{r) = \r, where A is a constant. The self-similar collapse 
solution is then 

R = (^-Aj r(l - zf/\ (3-94) 
and the partial derivative R with respect to r is 

where z = t/r. 

In the following, we derive the mapping function for the self-similar dust collapse 
acceding to Barve et al. ^^^^ For self-similar dust collapse, we can obtain exact null 
coordinates. For this purpose, we restrict our attention to the two-dimensional part 
of metric tensor. In the interior region, we introduce null coordinates (77, Q 

I „/"'/('-".■■) (z-R,>0) 



and 



We introduce another null coordinates {U, V) as 

^ = , f ^ nl (3-98) 

I - In \rj\ {z -R,r<0) ^ ^ 

\ -InlCI {z + R,r<0). ^^^^^ 

Hereafter we concentrate our attention to the outside of the Cauchy horizon where 
z — R^r < 0. Then the two-dimensional metric is expressed as 

ds'^ = -A^{U, V)dUdV (3-100) 

where 

If we define 

f±iy)^y'T^y'-yTla (3-102) 
1/3 



where y = {1 — z)^/^ and a=(^^Xj , we can write 

3r 
f±{y) 



r dz f Sy^dy 
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In the exterior region, the metric can be expressed as 



1 



— — ^ dudv 
R J 



(3-104) 



using the Eddington-Finkelstein double null coordinates. The matching surface is 
given by r = rj,. Hence, we have 2M = Ar;,. As in the analytic case, the matching 
condition is given by Eq. (1-31). 

3.3.2. Mapping 

The mapping v = G{u) from to I~ can be determined as follows. First we 
consider the limit to the regular center t < and r = 0. Since y ^ oo in this limit, 
we have 



T] -ry 



t, 
t. 



(3-105) 
(3-106) 



Therefore, rj = ( holds at the regular center. The singularity at t = and r = is 
mapped into rj = ( = 0. 

Next we consider the earliest outgoing null geodesic n_ which emanates from 
the center at i = 0. Since r] is outgoing null coordinate, rj = along n_. Since it 
is only possible when 7_ = — oo along the null geodesic, we find that y along n_ 
satisfies an algebraic equation f-{y) = 0. In the domain < y < 1, there exists two 
real roots for A < 6(26 — 15\/3), one degenerate root for A = 6(26 — 15\/3), and no 
real roots otherwise. A < 6(26 — 15\/3) is the condition for the naked singularity 
to exist in this model. It is clear that the largest positive root y = y_ G (0, 1) 
corresponds to n_, i.e., the Cauchy horizon. Hereafter we restrict our attention to 
the nondcgenerate case. For the first ingoing null ray n+ which terminates at the 
singular center at i = 0, we encounter a similar situation. In the domain y > 1, 
an algebraic equation /+(y) = has always the only one real root y = y+, which 
corresponds to rij^. 

Then we consider the matching condition for null rays. The matching condition 
at the surface is given in terms of u and v as 



/ ^ n 3 ^ 282, 

u = ub{y) = -—y - -^any - ny - rb in 



vb{y) 



a 

n 3 
—y 

a 



4 
3' 

4 282, 
-ar^y + r^y + -a In 



3y 
2a 
3y 
2a 



- 1 

- 1 



(3-107) 

(3-108) 



It is noted that these functions are regular around y = y±. 

We shall consider a family of null rays which emanate from X~ , cross the regular 
center, and reach X+. For an ingoing null ray which is earlier than and close to n+, 
from Eqs. (3-97) and (3-103), we can find 

C~-Kym-y+)^+, (3-109) 

where yj„ is the value of y with which the ingoing null ray crosses the surface and 

3yi 



7+ = 



f'+iy+y 



(3-110) 
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Similarly, for an outgoing null ray which is earlier than and close to n_, from 
Eqs. (3-96) and (3-103), we can find 

V ^ -riVout - y-y~ , (3-111) 
where yout is the value of y with which the ingoing null ray crosses the surface and 



7- = 



f'-{y- 



For an ingoing null ray which is earlier than and close to n+, we find 



V ^ Vq + 



dvb \ 
dy) 



{Vin -y+), 



(3-112) 



(3-113) 



y+ 



where vq = Vb{y+). Similarly, for an outgoing null ray which is earlier than and close 
to ri_, we find 



dub \ 

dy )y_ 



{yout -y-), 



(3-114) 



where uq = Ub{y-). The coefficients {dub/dy)y_ and {dvb/dy)y_^ are both finite and 
negative. Then we can find the following mapping from I'^ to I" using the fact that 
the regular center is given hy tj = (: 



V = G{u) vq + r{uo — uy, 



where 



7 = 



7- 



7+ 



r = 



dvb \ 
dy) 



a+ 



I dub \ 
\dy)o^_ 



-7 



(3-115) 

(3-116) 
(3-117) 



It should be noted that we can show 7 > 1. 

3.3.3. Power and energy 

Now that we have the map G for late times, we can calculate the emitted power 
and energy for late times. We have to assume that geometrical optics approximation 
is valid for a four-dimensional problem. The result is 



and 



P -- 
P-- 

E 
E = 



1 



7 



487r {uo - uf ' 

1 (7^1)^ 
487r {uq - uY ' 



1 1' 



1 



487r uq — u ' 

1 (7-1)^ ^ 
487r uq — u 



(3-118) 
(3-119) 

(3-120) 
(3-121) 
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Therefore, the emitted power diverges to positive infinity for both minimally and 
conformally coupled scalar fields. The totally radiated energy also diverges to pos- 
itive infinity for both minimally and conformally coupled scalar fields. It should be 
noted that there appears no characteristic scale in the above expression for emitted 
power and energy in contrast to the analytic case. 

We should note that no complete analysis for the spectrum of radiation has 
been done yet for the particle creation from self-similar dust collapse, although it 
was pointed out that the spectrum is different from that of a black body. ^^^^ 

3.3.4. Quantum stress-energy tensor 

Since we have an exact expression for the metric tensor in double null coordi- 
nates, we can estimate the vacuum expectation value of stress-energy tensor in a 
two-dimensional spacetime. ^^^^ 

The expectation value of stress-energy tensor in the interior coordinates is given 
by Eqs. (3-30)-(3-32). It is found that {Tjjy) in the interior is given by 



1 - Rl\d'^R,r 

r 

After a long calculation, we obtain 



127r 4 



127r4(l -Fx) 



where x is determined by the equation 




A 3x3 A 3^2y 



(3-123) 



^ = A + --A-21n^- (3-124) 
2M 3x3 a; x2 ^ ^ ' 

Fu{[3') can be expressed by the same equation (3-123) but with x which is determined 
by the equation related to the retarded time as 

r,-/3(a(.))^ 2 2_ 1 1 + x 

2M 3x3 X x2 X ^ ' 

Fu{A^) and Fv{A^) become 

^/.2x 1 (ifdRr V Iz'^ - Rld'^Rr] 
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In the exterior Schwarzschild region, the vacuum expectation value is given by 
Eqs. (3-35)~(3-37). We obtain 



1 

12^ 



X^-2A + 3 



+ — ^ fSw^ - 7w^) 

16 (2M)^ V ; 



127r 16(2M)2(3u;3 - 2A - 3u;4 - wX)^ 



a' (w'^ 2w^ 2w^ 
^ 4(2M) I T ^ Y 



X 3u;^A 



-(5^^ + 12u;^ - Sw^ - 2Aw + 12) A - 6w;^ + 18w^ - 12w;^l(3-128) 



where w is determined by the equation 

rb — u 2 



2 1 1-w 
+ - + ^ + 21n- 



2M 3w^ w w 

The relation between u and U is obtained as 

1 _ 2M /I 2 w 
a'{u) 1 — w \A 3w^ X 



w 



3w^ 



(3-129) 



(3-130) 



Here we interpret the outgoing flux emitted from the star. For this sake, we look 
into the right hand side of Eq. (3-35). The first term clearly denotes the vacuum 
polarization, which tends to vanish so rapidly as R goes to infinity that there is no 
net contribution to the flux at inflnity. The second term is originated in the interior 
of the star as seen in Eqs. (3-30) and (3-31). First it appears as ingoing flux at 
the passage of the ingoing rays through the stellar surface, crosses the center and 
becomes outgoing flux. Since this term depends on both a' and /?', not only the 
outgoing but also ingoing null rays are relevant to this term. It implies that this 
term strongly depends on the details of the spacetime geometry in the interior of the 
star. Since the third term depends only on a', only the outgoing null rays determine 
this contribution. Such a term is not seen in Eq. (3-30). Therefore, this contribution 
seems to originate at the passage of the outgoing rays through the stellar surface. As 
will be discussed later, the third term corresponds to the Hawking radiation, while 
the second term becomes important in the naked singularity explosion. 

3.3.5. Physical mechanism of particle creation 

Using the above exact expressions, we can now investigate the exact behaviors 
of the quantum stress-energy tensor for a massless scalar field. 

To determine what would be actually measured, the world line of an observer 
must be specified. For the observer with the velocity u^, the energy density {T^ 



and energy flux [T^i^u'^'ny are considered to be measured, where u'^n^ = and 



n^^n^ = 1. 



To investigate the importance of the back reaction for the central singularity 
formation, we compare the energy density observed by the comoving observer to the 
background energy density around the center. The energy density observed by the 
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(a) Lines of pq = const (b) Lines of p = const and pq = p 

Fig. 16. Plots of (a) lines pq = const and (b) p = const and pq = p in the self-similar LTB spacetime. 
Solid, long-dashed-dotted, dashed, dashed-dotted, long-dashed, long-dashed-short-dashed and 
dotted lines denote singularity, Cauchy horizon, pq = p(— 1,0), pq = 0, pq = — p(— 1,0), p = 
p(— 1,0) and pq = p, respectively. 



comoving observer becomes 

o ^(T )u''u^ - ^^^^^ I ^^^^^ t2 ^^^^^ f3-13n 

The results are shown in Fig. 16. Basically we use the background parameters A = 0.1 
and Tft = 10^. The lines of pq = const are plotted in (a) and the line of r) = const 
is in (b). These two values coincide with each other on the dotted line in (b). Below 
this line the background density is larger than the quantum energy density. We 
see especially that near the center the background density is larger than the energy 
density of the scalar field until the central singularity. Therefore we conclude that 
the back reaction does not become significant during the semi-classical evolution. 

Next we consider the energy flux measured by the comoving observer. For this 
observer, the energy flux becomes 

In the interior region, outgoing part of flux is proportional to {Tuu) and ingoing 
part of flux is proportional to (Tyy). Here we define the ingoing part Fin and the 
outgoing part Fout of the flux as 

^ {Tyy) ^ {Tuu) i qq^ 
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200 400 600 800 1000 
r 

Fig. 17. Signs of ingoing and outgoing parts of flux in the self-similar LTB spacetimc. The solid and 
long-dashed-dotted lines correspond to the singularity and the Cauchy horizon, respectively. The 
dashed line denotes the line on which the ingoing part of flux vanishes. On the long dashed lines 
the outgoing part vanishes. Regions 1, 11, 111, and IV correspond to positive ingoing and outgoing 
parts, positive ingoing and negative outgoing parts, negative ingoing and outgoing parts, and 
negative ingoing and positive outgoing parts, respectively. On the short-dashed-dotted line, the 
observed net flux vanishes. 



Following the sign of these two values we divide the interior region into four parts. 
Region I: positive ingoing and outgoing parts, region II: negative outgoing and pos- 
itive ingoing parts, region III: negative ingoing and outgoing parts, and region IV: 
positive outgoing and negative ingoing parts. The results are shown in Fig. 17. Flux 
observed by a time-like observer may become zero in region I and III. In fact, the 
observed flux vanishes at r = and on the dotted dashed line in Fig. 17. 

We show this flux schematically in Fig. 18(a). The intensity is proportional to 
the length of arrows along the r-axis. The comoving observers receive inward flux 
first and then outward flux after the crossing of the line (c). The intensity grows as 
the Cauchy horizon is approached. 

In Fig. 18(b), we show ingoing and outgoing parts of the flux schematically. The 
intensity is proportional to the length of arrow along the r-axis. The future and 
past directed arrows correspond to the positive and negative flux respectively. The 
outgoing part has large intensity near the Cauchy horizon. The ingoing part has 
large intensity only near the central singularity. 

Here we summarize the behavior of the quantum field in the interior of the star. 
At first inward flux appears and then the positive energy is concentrated near the 
center surrounded by a slightly negative energy region. As the collapse proceeds 
the central positive energy increases and concentrates within a smaller region. At 
the naked singularity formation the concentrated positive energy is converted to 
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(a) Observed flux (b) Ingoing and outgoing parts of flux 

Fig. 18. Schematic figures (a) for the flux observed by comoving observer and (b) for the ingoing 
and outgoing parts of flux in the self-similar LTB spacetime. Lines (a) and (b) correspond to 
the singularity and the Cauchy horizon, respectively. Line (c) denotes the line on which the 
observed flux vanishes. Regions I, II, III, and IV are similar to those in Fig. 17. 



the outgoing positive diverging flux. The left negative energy goes down across the 
Cauchy horizon. It will be reach the space-like singularity. 

Next we consider the totally radiated energy received at infinity. The radiated 
power of quantum flux can be expressed as 

P = {T^u)u^n'' = a'^FuiP') + F„(a') (3-134) 

where u'^ is considered as the velocity of a static observer. We plot this in Fig. 19(a). 
We reproduce the results of the geometric optics estimates for the four-dimensional 
model, i.e., the inverse square dependence for the retarded time in an approach to 
the Cauchy horizon. 

The totally radiated energy is obtained from the integration of P with respect 
to the retarded time at the future null infinity. The results are shown in Fig. 19(b). 
We see that the totally emitted energy is much smaller than the mass of the original 
star in the range that the semi-classical approximation can be trusted. We can see 
that the positive divergence of the radiated power comes from the term a''^Ftj(P') 
which has strong dependence on the spacetime geometry in the interior of the star. 

Let us turn to the inspection for the black hole formation. For larger u there 
appears difference between the naked singularity case and the covered one. This 
difference mainly depend on the behavior of a'. In an approach to the event horizon, 
w goes to unity and then a' goes to zero. It is easy to see that only the second term 



64 



T. Harada, H. Iguchi and K. Nakao 



10 







-10 



-20 



11 1 1 1 1 1 1 1 1 1 

sN\x=o 1 


1 1 1 1 1 1 1 1 1 1 1 




A=io-3 ; 


X=10-5\O^ 






, . . 1 . . . 1 w 


. . , 1 . . . 1 . . . 



-2 2 
log(u„-u) 

(a) Power 








-4-2 2 
log(u„-u) 

(b) Energy 



Fig. 19. (a) Radiated power and (b) integrated energy received at infinity in the self-similar LTB 
spacetime. Upper, middle and lower lines correspond to A = 0.1 ,10~^ and 10~^, respectively. 
Solid, dashed and long-dashed lines correspond to ri, = 10, 10^ and 10^, respectively. 



of the radiated power is left in this hmit and it becomes 

F„(a') = ^. (3-135) 

This means that F^^a') is interpreted as the Hawking radiation contribution. Prom 
Eq. (3-39) we obtain 

{Tvv) = o (3-136) 

^ ' 1927r(2M)2 ^ ' 

at the event horizon. This ingoing negative flux crossing into the black hole is 
balanced with the energy loss by the positive Hawking flux. 

There is similar balance between ingoing and outgoing parts at the Cauchy hori- 
zon for the naked singularity explosion. In the asymptotic region positive diverging 
flux is observed in an approach to the Cauchy horizon. It has been shown above 
that there is negative ingoing flux crossing the Cauchy horizon. This negative flux 
diverge only at the central naJied singularity. This diverging negative ingoing flux 
is balanced with the diverging positive outgoing flux which propagates along the 
Cauchy horizon. 

3.4. Naked singularities and quantum gravity 

The above discussions are based on quantum field theory in curved spacetime. It 
may lose its validity if we consider a time interval which is shorter than the Planck 
time t Planck- Therefore, the expressions for the energy radiated from a forming 
naked singularity can be conceivable at least at u ^ uq — tpianck- For self-similar 
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dust collapse, it is easily found that the emitted energy does not exceed the Planck 
energy ~ lO^^erg until this moment. For the analytic dust collapse, the discussion 
is not so straightforward, but the result is the same. ^^^^ 

It is noted that ujg may exceed the Planck frequency uJpianck because Ug is de- 
termined by some combination of two macroscopic scales, the free-fall time and the 
scale of inhomogeneity. Nevertheless, we can find that the totally radiated energy 
for this late-time radiation cannot exceed the Planck energy before the Planck time. 
The proof is the following. If uJs > ^pianck-, then it is found that the late-time ra- 
diation is only for uq — u < uj~^ < tpianck- Therefore, we only have to consider 
a;^ < uj Planck- this case, it is found 

^uo-u^tpianck ~ s'-^ PlanckY ^ < ^Planck, (3-137) 

^uo~u-^tpianck ^ '^Planck i In ) < (^Planck, (3-138) 

\^ Planck J 

where the second inequality comes from the fact that the function — xlnx has a 
maximum value in the domain < x < 1. Therefore, we can conclude that the 
totally radiated energy of the late-time emission does not exceed the Planck energy. 

It follows that the further evolution of the star will depend on quantum gravita- 
tional effects, and without invoking quantum gravity it is not possible to say whether 
the star radiates away on a short time or settles down into a black hole state. 

3.5. Discussion 

Quantum particle creation during the formation process of naked singularity 
in dust collapse results in diverging flux emission within a semi-classical theory of 
quantum fields in curved space. The discussions so far have been based on some 
specific models of naked singularities. However, calculations are dependent not on 
the matter fields but on the spacetime geometry around the singularity. Therefore, 
we can conjecture that the discussions and results will apply to a rather wider class 
of naked singularities which develop from regular initial data. Although it is not 
certain whether this phenomenon may be relevant to any observations, several au- 
thors have speculated naked singularities as a possible central engine of gamma-ray 
bursts. Any way, to answer the question may be out of scope of a semi- 

classical theory of quantum fields in curved space. Apart from the observational 
aspects, the contribution of emitted quantum particles as a source of gravitational 
field may become important in more realistic situations. The argument that it may 
prohibit the formation of naked singularities leads to the quantum version of cos- 
mic censorship conjecture. In particular, the balance relation between positively 
diverging outgoing flux to infinity and negatively diverging ingoing flux to naked 
singularity may suggest the quantum version of "naked singularity evaporation con- 
jecture" proposed by Nakamura, Shibata and Nakao. 

§4. Summary 

We have several examples in which naked singularity forms as a result of time 
development of regular initial Cauchy data, as we have seen in §1. However, we 
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should not conclude that they alone disprove the cosmic censorship conjecture in 
the exact meaning. In all of these examples, some kind of high symmetries, such 
as spherical symmetry, axisymmetry, cylindrical symmetry and self-similarity, arc 
assumed. Our knowledge on the effects of deviations from these symmetries on the 
naked singularities has been quite restricted so far. In critical collapse, it is known 
that the occurrence of a zero-mass black hole, which is regarded as naked singularity, 
is realized as a result of the exact fine-tuning. This means that it is a zero-measure 
event even within the assumed symmetry. In other cases, it is quite difficult in general 
to answer the question whether the appearance of naked singularity is generic or not. 
It should be also noted that there is no generally accepted standard measure of the 
space of initial data sets. Apart from genericity, in some of the examples of naked 
singularities, matter fields cannot be considered as realistic, such as a dust fluid and 
null dust. It is also difficult to answer the question what matter model is physically 
reasonable. Some people consider that elementary fields such as a scalar field will 
be physically reasonable matter models. However, generally speaking, it is only at 
the stage in which quantum effects become so important that such field descriptions 
of matter fields will be useful. Others consider that a system of infinitely many self- 
gravitating collisionless particles, i.e., the Einstein- Vlasov system will be suitable for 
the proof or disproof of the cosmic censorship conjecture. This proposition seems to 
be reasonable from the standing point that singularities that form also in Newtonian 
gravity should not be taken seriously, because of the global existence theorem of 
regular solutions with nonsingular initial data sets in the corresponding Newtonian 
system. However, unfortunately, we do not know any realistic matter field which 
is expected to behave like the Einstein- Vlasov system in an extremely high-density 
regime. 

However, these examples of naked singularities strongly motivate us to imagine, 
as a possible stage of gravitational collapse, extremely high-density and/or high- 
curvature regions uncovered by horizons. In these strongly curved regions, general 
relativistic effects will become crucial. In fact, it is considered that classical general 
relativity cannot apply beyond some energy scale, which will be characterized by 
the Planck scale. Therefore, the existence of classical solutions with naked singu- 
larity formed from regular initial data will imply not only the limitation of future 
predictability of general relativity but also the possibility that new frontier beyond 
general relativity may be exposed to us. Of course, if a naked singularity is a generic 
outcome of gravitational collapse, the importance to examine its consequences is 
much more increased. 

If naked singularity is a possible outcome of gravitational collapse, it may provide 
an exotic source of gravitational wave burst. Gravitational wave generation reflects 
the dynamical motion of the spacetime in strongly curved regions and they will 
directly reach us without obstructions unlike electromagnetic waves. As we have 
seen in §2, for linearized nonspherical perturbations of spherical dust collapse, not the 
energy flux of gravitational waves but the perturbation of the spacetime curvature 
blows up unboundedly as one approaches the Cauchy horizon of the background 
spacetime. This conclusion is also derived semi-analytically using the quadrupole 
formula for the gravitational wave emission from the dust collapse in Newtonian 
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gravity. From this fact, we may imagine that some sort of perturbations concerning 
spacetime curvature will be considerably amplified near the Cauchy horizon of a 
forming naked singularity, although it is not clear if this amplification can be related 
to any observable of a distant observer. This result also implies weak instability of 
the Cauchy horizon. Since it means the breakdown of the linearized perturbation 
scheme, some nonlinear analysis is necessary to make the conclusion clear. As we are 
having high sensitivity to gravitational waves by laser interferometric gravitational 
wave detectors, gravitational wave observations will provide interesting information 
of strongly curved regions. 

In addition to classical gravitational radiation, we have seen quantum particle 
creation from a forming naked singularity in the spherical dust collapse in §3. The 
calculation of this emission is based on quantTim field theory in curved spacetime. 
The calculation goes parallel to that of the Hawking radiation of black holes. The 
only difference is the form of the mapping function between the advanced and re- 
tarded time coordinates, which comes from the behavior of ingoing and outgoing 
radial null geodesies in the curved spacetime. The result is that a forming naked 
singularity emits explosive radiation at the final epoch. This conclusion comes from 
the fact that the mapping function has only lower differentiability at the Cauchy 
horizon. It is expected that this property of the mapping function will be common 
to a wide class of naked singularities. In the case of the two dimensional self-similar 
dust collapse, the back reaction effects of quantum radiation on the spacetime geom- 
etry will be negligible at least until the singularity formation. The radiated positive 
energy originates from the vacuum energy in the region surrounding the center. A 
large amount of positive energy is gradually concentrated and accumulated around 
the center in the course of the singularity formation. After that, the positive energy 
turns into positive outgoing flux just before the singularity formation. In addition 
to the back reaction effects, there appears a rather subtle problem between quantum 
field theory and quantum gravity. If we accept the Planck length a priori as the 
shortest length scale of quantum field theory, the total amount of radiated energy 
from a forming naked singularity cannot exceed the Planck energy as long as quan- 
tum field theory is justified. Actually, the total amount of radiated energy is greatly 
dependent on the shortest length scale of quantum field theory. 

In conclusion, the cosmic censorship conjecture has provided a strong motivation 
for researchers in this field. This conjecture was intended to conserve the complete- 
ness of future predictability of classical theory. However, we are now convinced of 
the limitation of classical theory as theory of everything. In this point of view, 
known examples of naked singularities imply that the future predictability of clas- 
sical theory is violated in a finite-measure set of initial data sets. Then, we may 
pay our attention to more general features of gravitational collapse. In fact, classical 
gravitational wave radiation and quantum particle emission from a forming naked 
singularity, which have been discussed in this paper, are both irrespective of the final 
outcome of the forming naked singularity, because these calculations are basically 
independent of the spacetime geometry in the causal future of the singularity. In 
this respect, the forming naked singularity can be regarded as an example of inho- 
mogeneously runaway collapsing systems. It could be said that our understanding of 
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gravitational collapse might have been too much based on the spherically symmetric 
homogeneous collapse, such as the Oppenheimer-Snyder model. Gravitational col- 
lapse in more general situations is now and will be an active and attractive field in 
gravitational physics. 
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Appendix A 

Gauge-Invariant Perturbations 

In this appendix we give a brief introduction to the formalism of Gerlach and 
Sengupta^^^' for perturbations around the most general spherically symmetric 
spacetime. 

We consider a general spherically symmetric spacetime with metric 

g^^dxf'dx" = gabix'^)dx''dx^ + R^ix'^)-fAB{x^)dx'^dx^ (A-1) 
and stress-energy tensor 

tf.^dx'^dx'' = tab{x'^)dx''dx'' + ^t/R^{x'^)jAB{x^)dx^dx^ , (A-2) 

where ^Asdx^dx^ = d9^ + sin"^ OdcfP. Lowercase Latin indices refer to radial and 
time coordinates, while uppercase Latin indices refer to 9 and (j). 

Now we introduce arbitrary perturbations of this spacetime. The angular de- 
pendence of perturbations is decomposed into series of tensorial spherical harmonics. 
The scalar spherical harmonics are Yi^{x^). A basis of vector harmonics is formed 
by Yi"!^ and 5™^ = (a^i^b- ^ basis of symmetric rank- two tensor harmonics is 
formed by Y^lAB. ^^is = ^/7AS + (^(^ + l)/2)^r7AB and 5n4:s+5f^B:.i- For i = 
and 1, the last two tensors vanish identically. Linear perturbations with different / 
and m decouple. In the following, we suppress these indices. We also suppress the 
explicit summation over them. Perturbations with different values of m for the same 
I have the same dynamics on a spherically symmetric background, and therefore m 
does not appear in the field equations. Spherical harmonics are called even parity 
if they have parity (—1)' under spatial inversion and odd parity if they have parity 
(— 1)'^^. Even and odd perturbations decouple. 

The odd-parity perturbations arc expressed as 

hf.^dxi'dx'' = ha{x'')SB{dx''dx^ + dx^dx") + h{x'')S^A:B)dx'^dx^ (A-3) 
for the metric and 

At^.^dxi'dx'' = Ata{x'')SA{dx''dx'^ + dx^dx") + At{x'')S^A:B)dx^dx^ (A-4) 
for matter. The even-parity perturbations are 

V = Kh{x'^)Ydx''dx^ + ha{x'^)Y.,B{dx''dx^ + dx^dx") 
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+[K{x'^)R^jabY + G{x'^)R^ZAB]dx^dx^ (A-5) 



for the metric and 



Atf,u = Atab{x'^)Ydx''dx'' + Ata{x'^)Y:B{dx''dx'' + dx^'dx") 

+[At^{x'^)R'^jABY + At'^{x'^)R^ZAB]dx^dx^ (A-6) 

for matter. Here covariant derivatives are distinguished as follows: 

JAB:C = 0, gab\c = 0- (A- 7) 

For convenience, we introduce 



and 



Va ^ ^, (A-8) 



Pa = ha- ^i?'G a. (A-9) 



We then introduce gauge-invariant variables to eliminate gauge ambiguities in 
perturbations. The gauge transformation is induced by the infinitesimal vector fields 

^^dx^ = M{x^)SAdx^, (A- 10) 

for odd parity and 

^^dx^" = e„(x")ydx" + ^{x^)YAdx^ (A-11) 
for even parity. The odd-parity gauge-invariant metric variables are given by 



ka = ha- -R'da l^j^ j . (A- 12) 

The odd-parity gauge-invariant matter variables are given by the combinations 

L^ = ta-\fj,^ha, (A- 13) 

L^t-\fj,^h. (A-14) 
A set of even-parity gauge-invariant metric perturbations is defined as 

Kb = Kb - {Pa\b + Pb\a), (A- 15) 

k = K+ i^ii-lic - 2vya- (A-16) 

A set of even-parity gauge-invariant matter perturbations is defined as 

Tab = ^tab - ta^cP" - taPc\b " h'^Pcla, (A- 17) 

Ta = Ata - t^pc - R\t//A)G,a, (A- 18) 

= At' - {pyR'){RH//2),, + l{l + m//A)G, (A-19) 

= At^ - {rH//2)G. (A-20) 
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The linearized Einstein equations, expressed only in terms of gauge-invariant 
perturbations, are 



My 

for odd parity and 



ka_^ _ j^A ( fee ^ 



= IGttL, (/ > 2) 



(A-21) 

{l-l){l + 2)ka = lQTTB?La, [l > X}.-22) 



(ykab\c kca\b fec6|a) 



i?2 



+ G,^ + G/ + 27^ 



+ i (G/ + G/) + 7^] fc/ + 2 (vak^ + Vbka + 



+9ab 



1(1 + 1) , 1 



i?2 



-dab 



2kJ'' + 6c''k 



(/-!)(/ + 2) 



i22 

fe,a ~ fe'ac' + |a ~ ''^a^c ~ "IGTrT^, 



167rr„6,(A-23) 
(A-24) 



-167rr3,(A-25) 



-167rr2 



(A-26) 

for even parity, where TZ is the Gaussian curvature of the 2-dimensional submanifold 
spanned by x", and 



Gab = -2 + VaV^ + 5„6 (^2v„l" + ^VaV"" - 

G/ = 2 (z;J'^ + z;,^;'^ - 7^) . 
The linearized conservation equation [6{T^^^) = 0] reduces to 
(i?2L«)^^ = (Z-l)(Z + 2)L, (Z>1) 

for odd parity and 
(i?2T»), 



(A-27) 
(A-28) 

(A-29) 



i?2 2i?2 

{B?Tabf TJil+1) 



T' = \t/ (fe - ^fe/) + ^f^Xb, (A-30) 



1 



i?2 



i?2 



- 2VaT' = h,,^at"^ + k^a - ^k^^^,t\ 



+2v%,t\ + k\t\^. 



(A-31) 



for even parity. 
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Appendix B 

Power of Gravitational Radiation 

In this appendix we examine the asymptotic behavior of the gauge-invariant 
variables in asymptotically flat spacetime with an outgoing wave condition. Then 
we can calculate the radiated power of the gravitational waves and thereby grasp 
the physical meaning of the gauge-invariant quantities. 

Note that in vacuum at large distance, the spherically symmetric background 
metric is identical to the Schwarzschild solution, which is the following (adopting 
the Schwarzschild coordinates): 

ds^ = _ (^1 _ dr^ + (^1 _ dR^ + (^^2 ^ ^^2 ^^^2 j _ (B.-L) 

To relate the perturbation of the metric to the radiated gravitational power, it is 
useful to specialize the gauge to the radiation gauge, in which the tetrad components 
^{0){e) ~ ^(<f)(<f) fs-ll off as 0(l/i?), and all other tetrad components fall 

off as 0{1/ B?) or faster with respect to the following background tetrad basis: 

= V~^) ^^'^^ 

ef,) = R{der, (B.4) 
el^^=Rsme{d^f. (B-5) 

In this radiation gauge, the metric perturbations in Eqs. (A-3) and (A-5) behave 

hQM = o(]-\, (B-6) 



for odd parity and as 



as 



R. 

h2 = w{t - R^) + 0{1) (B-7) 



hab = O (^-1) , (B-8) 
ha = (i) , (B-9) 
K = o(^), (B.IO) 



for even parity, where 



i?* = i?-F2M In - 1^ const, (B-12) 
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and the outgoing wave condition is satisfied. Then, the gauge-invariant metric per- 
turbations (A- 15) and (A- 16) are calculated as 



ko 



lw^^'^R + 0{l), (B-13) 



2 



ki = ^w^^^R + 0{l) (B-14) 

for odd parity and 

krT = g^^'^R + Oil), (B-15) 

krR = -9^^^R + 0{l), (B-16) 

kRR=g^^'>R + 0{l), (B-17) 
fc = -5« + 0(i) (B-18) 

for even parity, where w^'"^ and g^'"'> denotes the n-th derivatives of w and g, respec- 
tively. 

In this radiation gauge, the radiated power P per soUd angle is given by the 
formula derived by Landau and Lifshitz from their stress-energy pseudo-tensor: 

dP _ R^ 



(B-iyj 



dr J 4\ dr dr J 



For the axisymmetric mode, (i.e. m = 0), the above formula is reduced to 

(B-20) 

for odd parity and 



dP _ 


1 


df2 ~ 




dP 


1 


dn ~ 


64^ 



(B-21) 



for even parity, where 



^ ' Att \ (dcos6')2 J ^ ^ 

It is found that for the monopole and dipole modes, the radiated power exactly 
vanishes. By using gauge-invariant quantities and integrating over the 47r solid angle, 
the formula for the power of the gravitational radiation for Z > 2 is obtained as 



for odd parity and 



dP 1 

P = —Bik^ (B-26) 
647r ^ ^ 
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for even parity, where 



R - (^ + ^)' fR 27^ 



Appendix C 

Equations for Even-Parity Perturbations 

For marginally bound LTB spacetime linearized quadrupole Einstein equations 

are 

4 1 , R 

~r^''^[rr~r^)^ ~ [ r^r)''^rJ^'' 

I R RR'' R' \ , R , , „ r- / ^ N 



R'^ 



R' 

r' 



R'^R 
R 



-g + 4- 



R'^R, 



R 



k + R'^k 



R' R R' • 



0, 



(C-3) 



-2§jq -q + 2^k + 2k + ^koi - -^koi' = -167rpy2, (C-4) 



R' 

q + koi + -^koi = 0, (C-5) 



R^R" R\ I / „^ R^R'\ R? ,1 „9 . 

2-\q'-[2RR + i—^]q-—^q"-R\ 



R'3 R'J^ \ R' i?'2 

/„• R^R'\ ■ „„2;- ^(RR"R RR'\ , 
+4 + j k + 2i?2fc + 2 - j /coi 

-2|^A.oi' + 2 ^ - I Uoi - = 0, (C-6) 



-feoo + = 0. (C-7) 

Here we can use Eq. (C-7) to eliminate k\i in Eqs. (C-l)-(C-6). 
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Appendix D 

Characteristic Frequency of a Naked Singularity 

In §1.2.2, we set the initial data at t = and specified the solution. At an 
arbitrary time t < to, we can expand the mass function F(r) in terms of the cir- 
cumferential radius R for marginally bound collapse with an analytic initial density 
profile as 

F{r) = Fa ( '-^) + fJ '-^) R' + ■ ■ ■ . (D-l) 



to J \ to 

Although the quantities F{r), {to — t) and R have physical meanings, to itself has 
no physical meaning, because it explicitly depends on the choice of the origin of the 
time coordinate t. This implies that neither F^ nor F5, but, rather, the combination 
F^^/F^ characterizes the time evolution of the density profile in the neighborhood 
of the center, because the dependence on to vanishes only through this combination 
of F3 and F5. 

Observing that the free-fall time and the scale of the inhomogeneity at the time 
t = are 

to ^ ^F3-^/^ (D-2) 
k - (^) , (D-3) 

we can easily find that the following constant, which has the dimension of frequency, 
characterizes naked singularity formation in the LTB spacetime: 

Us ^ |. (D.4) 



References 

R. Penrose, Phys. Rev. Lett. 14 (1965), 57. 

S. W. Hawking, Proc. Roy. Soc. London A 300 (1967), 187. 

S. W. Hawking and R. Penrose, Proc. Roy. Soc. London A 314 (1970), 529. 

R. Penrose, Riv. Nuovo Cim. 1 (1969), 252. 

R. Penrose, in General Relativity, an Einstein Centenary Survey, ed. S.W. Hawking and 
W. Israel (Cambridge University Press, Cambridge, England, 1979), p. 581. 
R. Casadio and B. Harms, hep-th/0110255. 

S. W. Hawking and G.F.R. Ellis, The large scale structure of space-time (Cambridge Uni- 
versity Press, Cambridge, England, 1973). 
C. W. Misner and D.H. Sharp, Phys. Rev. 136 (1964), 571. 

F. J. Tipler, Phys. Lett. A 64 (1977), 8. 
A. Krolak, J. Math. Phys. 28 (1987), 138. 

C. J. S. Clarke, The Analysis of Space-Time Singularities (Cambridge University Press, 

Cambridge, England, 1993). 

G. Lcmaitrc, Ann. Soc. Sci. BruxcUcs 1 A 53 (1933), 51. 
R. C. Tolman, Proc. Nat. Acad. Sci. 20 (1934), 169. 

H. Bondi, Mon. Not. R. Astron. Soc. 107 (1947), 410. 

P. Yodzis, H.-J. Seifert and H. Miiller zum Hagen, Commun. Math. Phys. 34 (1973), 135. 



Physical Processes in Naked Singularity Formation 



75 



16) P. S. Joshi and I. H. Dwivedi, Phys. Rev. D 45 (1993), 5357. 

17) D. M. Eardley and L. Smarr, Phys. Rev. D 19 (1979), 2239. 

18) D. Christodoulou, Commun. Math. Phys. 93 (1984), 171. 

19) T. P. Singh and P. S. Joshi, Class. Quantum Grav. 13 (1996), 559. 

20) S. Jhingan and P. S. Joshi, Annals of Israel Physical Society Vol. 13 (1997), 357. 

21) R. P. A. C. Newman, Class. Quantum Grav. 3 (1986), 527. 

22) S. S. Doshingkar, P. S. Joshi and 1. H. Dwivedi, Phys. Rev. D 59 (1999), 044018. 

23) I. H. Dwivedi, Phys. Rev. D 58 (1998), 064004. 

24) F. Mena and B. Nolan, Class. Quantum Grav. 18 (2001), 4531. 

25) S. S. Deshingkar, P. S. Joshi and I. H. Dwivedi, gr-qc/0111053. 

26) J. R. Oppenheimer and H. Snyder, Phys. Rev. 56 (1939), 455. 

27) H. Miiller zum Hagen, P. Yodzis and H. -J. Seifert, Commun. Math. Phys. 37 (1974), 29. 

28) A. Ori and T. Piran, Phys. Rev. Lett. 59 (1987), 2137. 

29) A. Ori and T. Piran, Gen. Relat. Gravit. 20 (1988), 7. 

30) A. Ori and T. Piran, Phys. Rev. D 42 (1990), 1068. 

31) T. Foglizzo and R. N. Henriksen, Phys. Rev. D 48 (1993), 4645. 

32) B. Waugh and K. Lake, Phys. Rev. D 38 (1988), 1315. 

33) B. Waugh and K. Lake, Phys. Rev. D 40 (1989), 2137. 

34) T. Harada, Phys. Rev. D 58 (1998), 104015. 

35) T. Harada and H. Maeda, Phys. Rev. D 63 (2001), 084022. 

36) T. Harada, Class. Quantum Grav. 18 (2001), 4549. 

37) B. J. Carr, gr-qc/0003009. 

38) B. J. Carr and A. A. Coley, Class. Quantum Grav. 16 (1999), R31. 

39) M. W. Choptuik, Phys. Rev. Lett. 70 (1993), 9. 

40) A. M. Abrahams and C. R. Evans, Phys. Rev. Lett. 70 (1993), 2980. 

41) C. R. Evans and J. S. Coleman, Phys. Rev. Lett. 72 (1994), 1782. 

42) D. W. Neilsen and M. W. Choptuik, Class. Quantum Grav. 17 (2000), 733. 

43) D. W. Neilsen and M. W. Choptuik, Class. Quantum Grav. 17 (2000), 761. 

44) T. Koike, T. Hara and S. Adachi, Phys. Rev. Lett. 74 (1995), 5170. 

45) H. Maeda and T. Harada, Phys. Rev. D 64 (2001), 124024. 

46) C. Gundlach, Liv. Rev. Rel. 2 (1999), 4. 

47) D. Christodoulou, Conmnm. Math. Phys. 109 (1987), 613. 

48) D. Christodoulou, Commun. Pure Appl. Math. 44 (1991), 339. 

49) D. Christodoulou, Commun. Pure Appl. Math. 46 (1993), 1131. 

50) D. Christodoulou, Ann. of Math. 140 (1997), 607. 

51) D. Christodoulou, Ann. of Math. 149 (1999), 187. 

52) A. Einstein, Preuss. Akad. Wiss. Berlin Sitzber. (1915), 844. 

53) B. K. Datta, Gen. Relat. Gravit. 1 (1970), 19. 

54) H. Bondi, Gen. Relat. Gravit. 2 (1971), 321. 

55) A. B. Evans, Gen. Relat. Gravit. 8 (1976), 155. 

56) G. Magh, Class. Quantum Grav. 14 (1997), 1937. 

57) G. Magh, Class. Quantum Grav. 15 (1998), 3215. 

58) T. Harada, K. Nakao and H. Iguchi, Class. Quantum Grav. 16 (1999), 2785. 

59) T. Harada, H. Iguchi and K. Nakao, Phys. Rev. D 58 (1998), 041502. 

60) S. Jhingan and G. Magh, Phys. Rev. D 61 (2000), 124006. 

61) H. Kudoh, T. Harada and H. Iguchi, Phys. Rev. D 62 (2000), 104016. 

62) A. D. Rcndall, Liv. Rev. Rel. 3 (2000), 1. 

63) K. Pfaffelmoser, J. Diff. Eq. 95 (1992), 281. 

64) A. D. Rendall, in Approaches to Numerical Relativity, ed. R.A. d'Inverno (Cambridge 
University Press, Cambridge, England, 1992), p. 94. 

65) G. Rein and A. D. Rendall, Commun. Math. Phys. 150 (1992), 561. 

66) G. Rein, A. D. Rendall and J. Schaeffer, Commun. Math. Phys. 168 (1995), 467. 

67) P. C. Vaidya, Current Science 13 (1943), 183. 

68) P. C. Vaidya, Phys. Rev. 47 (1951), 10. 

69) P. C. Vaidya, Proc. Indian Acad. Sci. A 33 (1951), 264. 

70) P. S. Joshi and 1. H. Dwivedi, Commun. Math. Phys. 146 (1992), 333. 

71) P. Szekeres, Phys. Rev. D 12 (1975), 2941. 

72) P. S. Joshi and A. Krolak, Class. Quantum Grav. 13 (1996), 3069. 



76 



T. Harada, H. Iguchi and K. Nakao 



73) S. S. Dcshingkar, S. Jhingan and P. S. Joshi, Gen. Rclat. Grav. 30 (1998), 1477. 

74) S. Chandrasekhar and J. B. Hartle, Proc. R. Soc. London A 384 (1982), 301. 

75) E. Poisson and W. Israel, Phys. Rev. D 41 (1990), 1796. 

76) P. R. Brady and E. Poisson, Class. Quantum Grav. 9 (1992), 121. 

77) P. R. Brady, I. G. Moss and R. C. Myers, Phys. Rev. Lett. 80 (1998), 3432. 

78) S. Hod and T. Piran, Phys. Rev. Lett. 81 (1998), 1554. 

79) A. Ori, Phys. Rev. D 58 (1998), 084016. 

80) A. Ori, Phys. Rev. Lett. 83 (1999), 5423. 

81) K. S. Thorne, in Magic Without Magic, ed. J.R. Klauder (Preeman, San Francisco, USA, 
1972), p 231. 

82) T. Nakamura and H. Sato, Prog. Theor. Phys. 67 (1982), 1396. 

83) T. Nakamura, K. Maeda, S. Miyama and M. Sasaki, in Proceedings of the 2nd Marcel 
Grossmann Meeting on General Relativity, ed. R. RufRni (North-Holand, Amsterdam, 
Netherland, 1982), p. 675. 

84) S. L. Shapiro and S. A. Tcukolsky, Phys. Rev. Lett. 66 (1991), 994. 

85) S. L. Shapiro and S. A. Tcukolsky, Phys. Rev. D 45 (1992), 2006. 

86) R. M. Wald and V. Iyer , Phys. Rev. D 44 (1991), R3719. 

87) C. Barrabes, W. Israel and P. S. Letelier, Phys. Lett. A 160 (1991), 41. 

88) M. A. Pelath, K. P. Tod and R. M. Wald, Class. Quantum Grav. 15 (1998), 3917. 

89) H. Iguchi, K. Nakao and T. Harada, Phys. Rev. D 57 (1998), 7262. 

90) H. Iguchi, T. Harada and K. Nakao, Prog. Theor. Phys. 101 (1999), 1235. 

91) U. H. Gcrlach and U. K. Scngupta, Phys. Rev. D 19 (1979), 2268. 

92) U. H. Gcrlach and U. K. Scngupta, Phys. Rev. D 22 (1980), 1300. 

93) J. M. Bardeen and T. Piran, Phys. Rep. 96 (1983), 205. 

94) S. Chandrasekhar and S. Dctweiler, Proc. R. Soc. London 344 (1975), 441. 

95) C. T. Cunningham, R. H. Price and V. Moncrief, Astrophys. J. 224 (1978), 643. 

96) C. T. Cunningham, R. H. Price and V. Moncrief, Astrophys. J. 230 (1979), 870. 

97) C. T. Cunningham, R. H. Price and V. Moncrief, Astrophys. .1. 236 (1980), 674. 

98) H. Iguchi, T. Harada and K. Nakao, Prog. Theor. Phys. 103 (2000), 53. 

99) T. Nakamura, M. Sliibata and K. Nakao, Prog. Theor. Phys. 89 (1993), 821. 

100) K. Nakao, H. Iguchi and T. Harada, Phys. Rev. D 63 (2001), 084003. 

101) T. A. Apostolatos and K. S. Thorne, Phys. Rev. D 46 (1992), 2435. 

102) F. Echeverria, Phys. Rev. D 47 (1993), 2271. 

103) T. Chiba, Prog. Theor. Phys. 95 (1996), 321. 

104) S. W. Hawking, Commun. Math. Phys. 43 (1975), 199. 

105) L. H. Ford and L. Parker, Phys. Rev. D 17 (1978), 1485. 

106) W. A. Hiscock, L. G. Wilhams and D. M. Eardley, Phys. Rev. D 26 (1982), 751. 

107) N. D. Birrcl and P. C. W. Davics, Quantum Fields m Curved Space (Cambridbe University 
Press, Cambridge, England, 1982). 

108) P. C. Davics, S. A. Fulling and W. G. Unruh, Phys. Rev. D 13 (1976), 2720. 

109) T. Harada, H. Iguchi and K. Nakao, Phys. Rev. D 61 (2000), 101502(R). 

110) T. Harada, H. Iguchi and K. Nakao, Phys. Rev. D 62 (2000), 084037. 

111) T. Tanaka and T.P. Singh, Phys. Rev. D 63 (2001), 124021. 

112) S. Barvc, T. P. Singh, C. Vaz and L. Wittcn, Nucl. Phys. B 532 (1998), 361. 

113) C. Vaz and L. Wittcn, Phys. Lett. B 442 (1998), 90. 

114) S. Barvc, T. P. Singh, C. Vaz and L. Witten, Phys. Rev. D 58 (1998), 104018. 

115) H. Iguchi and T. Harada, Class. Quantum Grav. 18 (2001), 3681. 

116) T. Harada, H. Iguchi, K. Nakao, T. P. Singh, T. Tanaka, and C. Vaz, Phys. Rev. D 64 
(2001), 041501(R). 

117) P. S. Joshi, N. K. Dadhich and R. Maartens, Mod. Phys. Lett. A 15 (2000), 991. 

118) L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields (Pergamon, London, 
England, 1975). 



